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ABSTRACT 

Aims. The Aim of our study is to understand the variety of observed Lya line profiles and strengths in Lyman Break Galaxies (LBGs) 
and Lya- emitters (LAEs), the physical parameters governing them, and hence to derive constraints on the gas and dust content and 
stellar populations of these objects. 

Methods. Using our 3D Lya radiation transfer code including gas and dust (Verhamme et al. 2006), we fit 1 1 LBGs from the FORS 
Deep Field with redshifts between 2.8 and 5 observed by Tapken et al. (2007). A simple geometry of a spherically expanding shell of 
H i is adopted. 

Results. The variety of observed Lya profiles is successfully reproduced. Most objects show outflow velocities of V cx( , ~ 150-200 
km s _1 ; two objects are most likely quasi-static. The radial H i column density ranges from N m ~ 2 x 10 19 to 7 x 10 20 cm . Our 
Lya profile fits yield values of E(B — V) ~ 0.05-0.2 for the gas extinction. We find indications for a dust-to-gas ratio higher than the 
Galactic value, and for a substantial scatter. The escape fraction of Lya- photons is found to be determined primarily by the extinction, 
and a simple fit formula is proposed. In this case a measurement of EW(Lya) b s can yield E(B - V), if the intrinsic Lya equivalent 
width is known (or assumed). Intrinsic EW(Lya)j nt ~ 50-100 A are found for 8/11 objects, as expected for stellar populations forming 
constantly over long periods (^ 10-100 Myr). In three cases we found indications for younger populations. Our model results allow 
us also to understand observed correlations between EW(Lya) obs and other observables such as FWHM(Lya), E(B - V), SFR(UV) 
etc. 

We suggest that most observed trends of Lya, both in LBGs and LAEs, are driven by variations of N ai and the accompanying variation 
of the dust content. Ultimately, the main parameter responsible for these variations may be the galaxy mass. We also show that there is 
a clear overlap between LBGs and LAEs: at z ~ 3 approximately 20-25 % of the LBGs of Shapley et al. (2003) overlap with ~ 23 % 
of the LAEs of Gronwal et al. (2007). Radiation transfer and dust effects should also naturally explain the increase of the LAE/LBG 
ratio, and a higher percentage of LBGs with strong Lya emission with increasing redshift. 

Key words. Galaxies: starburst - Galaxies: ISM - Galaxies: high-redshift - Ultraviolet: galaxies - Radiative transfer - Line: profiles 
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^ . 1. Introduction ionisation of the intergalactic medium and hence on cosmic 

' reionisati on, trace l arge scale structure at high redshift etc. (see 

Lya line radiation, often of the brightest emission lines in e , g J S chaerdl2007i for an overview, and references therein). 



distant star-forming galaxies, is nowfrequently observedoyer Howeveri given the physics of this generally optically thick res- 
a wide redshift ranggjHu et al. jl928k JCudntzki et al. 2000; Qnance Hne quantitative interpretations of Lya are often difficult 
Malhotra & Rhoads 2002; Ajiki et al. 2003; Tamguchi et al. nr _ v „ atTlWllnlls sWn hv A Pfai UA st]1 H; P . s nf n^rhv s tar- 



Aji kietalJ |2003|; |Taniguchi et al] or eyen am bi gU ous, as shown by detailed s tudies of n earby star- 
2005; Shimasakuetal.1 |2006j pshikawa^talj \200A bursts carr ied out during the last de cade (iLeaueux et alJ ll995; 
Tapken etalj |2006k | Gronwall et alJ |2007j; l Ouchietal] |2007l r^ unth et all S iMas-Hesseltd] 120031 lHaves et alT 12005. 



e.g.). Numerous narrow-band and other surveys use this line to 20071; I Atek et a l 2008) 
search for galaxies at specific cosmic ages, and to hunt for the 
most distant objec t s in the universe Xsee lWillis & Courbinl2005t 

ICubv et alJ 120071; IStark et alJ l2007al) . Furthermore extremely With the increased computer po wer and the avai l ability 
deep "blind" spectroscopic exposures h ave revealed ver y of Lya radiation transfer codes (cf. lAhn et all l2000l |200: 



faint objects through their Lya emission dRauch et alJ 120071) . Cantalupo et al. 2005; Hansen & Oh 2006; Diikstra et al]2006at 

illustrating also the discovery potential of future instruments iTasitsiomil l2006t IVerhamme et a ll 120061) time is now ripe not 

such as the Multi Unit Spectroscopic Explorer (MUSE) for the only to predict Lya in different astrophysical situations (e.g. 

Very Large Telescope (VLT) and extremely large telescopes. galaxy and cosmological simulations) but especially also to con- 

Lya measurements are used to infer a number of proper- front observed Lya properties (line profiles, equivalent widths, 

ties such as redshift, star formation rates, constraints on the etc.) of individual galaxies (nearby and distant ones) with de- 

tailed 3D radiation transfer calculations. This is one of the aims 

Send offprint requests to: daniel.schaerer@obs.unige.ch of our paper. 
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Two galaxy populations, the well known Lyman Break 
Galaxies (LBGs) and the Lye emitters (LAEs), represent cur- 
rently the largest samples of distant galaxies, at least from z ~ 
3 to 6.5. Important questions remain unanswered about them, 
closely related to their Lya emission and absorption. 

LBGs show a great diversity of Lya profiles and strengths, 
reaching from strong emissi on, over P-Cygni ty pe profiles, to 
strong absorption (see e.g. IShaplev et al.l 120031) . Furthermore 
the Lya properties are found to correlate with other quanti- 
ties such as the strength of interstellar lines (IS), extinction, 
the star formation rate (SFR) etc. However, the origin of these 
variations and correlations re mains largely unknow n or contra- 
dictory at best. For example, IShaplev et al.l (12001 1) suggested, 
age as the main difference between LBG groups with differ- 
ent Lya strength, where objects with Lya in absorption would 
be younger and more dusty. However, it is not clear why older 
LBGs would contain less dust, especially since outflows, sup- 
posedly used to expel the dust, are ubiquitous in all LBGs 
and since precisely these outflows ar e the location where the 
emergent Lya spectrum is determined. iFerrara &Ricottil (120061) 
proposed that LBGs host short-lived (30 + 5 Myr) starburst 
episodes, whose outflows - when observed at different evolu- 
tionary phases - would give rise to the observed correlations be- 
tween IS lines and Lya. However, the ages obtained from SED 
and spectral fits of LBGs s how older ages and fairly constant 
star formation histories ( cf. [Ellingson et al.l 1 1996t iPe ttini et al. 
2000; IShaplev et ai1l200ll;lPapovich et alJl2001l;lPentericci et alj 



2007). lPentericci et all (12007) find that z ~ 4 LBGs without Lya 
emission are on average somewhat older (~ 410 ± 70 versus 
200 + 50 Myr) and more massive than those with Lya in emis- 
sion. However, given the significant amount of ongoing star for- 
mation inferred for both galaxy subsamples (with/without Lya 
emission) types from their SED fits, intrinsic Lya emission is 
expected in both subsamples. The apparent differences in age 
and ratio of age over star formation time scale (t) can therefore 
not be a physical cause for the observed Lya variations. lErb et all 
(2006), from an analysis of z ~ 2 LBGs, also find that more mas- 
sive galaxies ha ve fainter Lya (or Lya i n absorption); following 
the arguments of M as-Hesse et al.l d2003l) they suggest that this is 
mostly due to an increased velocity dispersion of the interstellar 
medium, indicated by the increased strength of the saturated IS 
lines, which would increase the fraction of Lya being absorbed. 
However, from radiation transfer modelling we find that the be- 
haviour of the Lya escap e fraction is not line ar with the velocity 
dispersion in the ISM (Verha mme et alj|2006l and below). Other 
physical parameters may play a more important role in the de- 
struction of Lya. 

Important questions remain also concerning the proper- 
ties of LAEs, their similarities and differences with respect 
to LBGs, and about the overlap between these galaxy pop- 
ulations. Since a fraction of LBGs show strong enough Lya 
emission to be detected with the narrow-band technique mostly 
used to find LAEs there must definitely be an overlap. For 
example, at z ~ 3 ap proximately 25 % of the LBGs of 
Steidel and collaborators ( Shapley et al. 20033 have EW(Lya) b s 
<: 20 A (restframe), sufficien t to be detected in the LAE sur- 
vey of iGronwall et al.l d2007j) . However, what the properties of 
LAEs are e.g. in terms of stellar populations (age, star for- 
mation histories, . . .), mass, dust content, outflows, metallic- 
ity etc. is not yet well established, alt hough first such analy- 
sis have recently become available (s ee Schaerer & Pello 2005; 
Lai et al.ll2007l 120081; IGronwall et alJl2007HGawiser et alj|2007t 
Pirzk al et a 1J I2007I; iFinkelstein et al.ll2007bl) Understanding the 
nature of LAEs and their relation to LBGs is also crucial since 



the contribution of the LAE population to t he known star- 
burst population seems t o increase with redshift dHu et al.l [l998; 
I Shim asaku et all 120061; iNagao et alj 120071; lOuchi et alj 120071: 
bow-Hygelund et alJl2007l;lReddv et al.ll2007l) . 

Last, but not least, several types of theoretical models 
have been constructed during the last few years aimed at 
understanding LAE and LBG populations, the relation be- 
tween the two populations, and to use them as constraints for 
galaxy formation scenarios, cosmi c reion is ation, and other 
topics dThommes & Meisenheimerj 120051; iLe Delliou et ail 
1 20061; iMori & Umemural 120061: piikstra & Wvithd 12007; 
Kobavashi et alj 120071; iMao et alj [20071; IStarketalJ [2007b; 
Nagamine et al. 2008). Observational constraints on crucial 
parameters such as the Lya escape fraction from LBGs and 
LAEs, and other insight from radiation transfer models are, 
however, badly needed to reduce uncertainties and degeneracies 
in these modeling approaches. 

With these questions about LBGs and LAEs in mind, we 
have recently started to model a variet y of z ~ 3 starbursts 
with our new Lya radiation transfer code ( Verha mme et alj2 006, 
paper I). First we have studied the well known z ~ 2.7 LBG 
MS1512-cB58 (cB58 in short), whose spectrum is dom inated by 
strong Lya absorption (Scha erer & Ver hamme 2008, paper II). 
In this third paper of the series we present an analysis of 11 
LBGs observed with FORS2 at the VLT with sufficient spec- 
tral resolution (R ~ 2000) to allow detailed Lya profile fitting to 
constrain their properties. Taken together, the objects analysed in 
paper II and III cover a wide range of Lya strengths and also dif- 
ferent morphologies, including absorption dominated Lya and 
Lya emission lines with equivalent widths between ~ 6 and 150 
A (restframe). The variety of objects modeled in paper II and III 
covers thus in particular the entire range of Lya s trengths defin- 
ing th e 4 spectral groups of the LBG sample of IShaplev et al.1 
d2003h . the largest currently available at z ~ 3. Furthermore, sev- 
eral of the objects we model have strong enough Lya emission to 
classify as LAEs, according to the criteria used in many surveys. 
Our analysis represents the first modeling attempt of Lya line 
profiles of high redshift galaxies with a detailed radiation trans- 
fer code including gas and dust and treating line and continuum 
radiation. 

The remainder of the paper is structured as follows. A de- 
scription of the radiation transfer code, the assumptions, and in- 
put parameters is given in Sect. [2] In Sect. [3] we model the Lya 
profiles of the individual objects. Our main fitting results are dis- 
cussed and confronted to other observations in Sect. [4] Other 
properties are derived in Sect. [5] In Sect. [6] we finally propose a 
unifying scenario for LBGs and LAEs and discuss several impli- 
cations. Our main conclusions are summarised in Sect. [7] 



2. Lyn radiation transfer modeling 

To fit the observations we us ed our 3D Monte Carl o (MC) ra- 
diation transfer code MCLya dVer hamme et al. 2006). The code 
solves the transfer of Lya line and adjacent continuum photons 
including the detailed processes of Lya line scattering, dust scat- 
tering, and dust absorption. The main assumptions required for 
the modeling concern the geometry, the choice of the input pa- 
rameters, and the input spectrum. We discuss them now in turn. 



2.1. Geometry 

For simplicity, and given empirical evidence in favour of a fairly 
simple geometry in z ~ 3 LBGs discussed by IVerhamme et al.l 
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(2006), we adopted first a simple "super-bubble" model to at- 
tempt to fit the observed Lya line profile. The assumed geometry 
is that of an expanding, spherical, homogeneous, and isothermal 
shell of neutral hydrogen surrounding a central starburst emitting 
a UV continuum plus Lya recombination line radiation from its 
associated H n region. We assume that dust and H i are uniformly 
mixed. 

The homogeneity and a large covering factor of the absorb- 
ing medium are supported by observations of LBGs and nearby 
starbursts. For example, the outflow of cB58 is located well in 
front of the stars and covers them almost completely, since it ab- 
sorbs almost a l l the U V light from the background stars. Indeed, 
ISavaglio et al. |(|2002|)find only a small residual mean flux above 
the zero level i n the core of t he Lya absorption featu r e, whe reas 
it is black for iPettini et all d2002l) . iHeckman etall (1200 ll) es- 
timate an area covering factor for optically thick gas of 98% 
from the residual intensity at the core of the C n ^11 335 line. 
A somewhat lower coveri ng factor may be ind icated for LBGs 
with strong Lya emission (Shapley et al. 2003). For simplicity, 
and in the absence of further observational constraints, we will 
assume a covering factor of unity. 

2.2. Shell parameters 

As described in paper II, the outflow is modelised by a spheri- 
cal, homogeneous, and isothermal shell of neutral hydrogen and 
dust centered on a point source. Four parameters characterise the 
physical conditions in the shell: 

- the expansion velocity Ve xp , 

- the Doppler parameter b, 

- the neutral hydrogen column density Nm, 

- the dust absorption optical depth r fl . 

In principle V e xp is constrained by observations, either di- 
rectly measured by the blueshift of low ionisatio n interstel- 
lar lin es (hereafter LIS) compared to stellar lines (Petti nTet alj 
|2002|) . or from the shift between absorption LIS lines and Lya 
in emission, Av(em - abs), when the stellar lines are too faint to 
be observed dShaplev et al.ll2003l) . Otherwise V exv w i ll be c on- 
strained by Lya line profile fits. In [Verhamme et al. (2006]) we 
showed that radiation transfer effects lead to Av(em - abs) « 
3 x V exp in expanding shells with Nui ^ 10 20 crrr 2 . For lower 
column densities the peak of the redshifted Lya emission may 
trace ~ V exp ( leading to Av(em — abs ) ~ 2V exp ), instead of twice 
this value (cf. Verhamme et al. 2006, and also FigfTTb. For three 
of the 1 1 objects to be modeled here Av(em- abs) has been mea- 
sured (see Table[TJ. 

The Doppler parameter b, describing the random motions of 
the neutral gaz possibly including microturbulence, is kept as 
a free parameter. For indication, b ~ 13 km s corresponds 
to thermal motion s for T = 10 4 K; for the lensed LBG cB58 
IPettini et al] d2002l) derived b ~ 70 km s _1 from fits of LIS lines. 

Although presumably the neutral column density and the 
dust amount are physically related, e.g. by a given dust-to-gas 
ratio, both parameters are kept free in our modeling procedure. 
The resulting values will later be compared to available observa- 
tional constraints. 

We as sume that dust and H 1 are uniformly mixed. As dis- 
cussed in IVerhamme et alj d2006l) . the dust optical depth j a re- 
lates to the usual extinction E(B - V) « (0.06 ... 0. 1 1) r a , where 
the numerica l coeffi c ient cov ers attenuation/extinction laws of 
ICalzetti et ail d2000t) . ISeatonl d 19791) and similar. Here we as- 
sume E(B-V) = 0.lT fl . 



2.3. The intrinsic spectrum in the Lya region 

The synthetic stellar spectrum of star-forming galaxies close to 
Lya is described in paper II. The stellar continuum presents an 
absorption feature around Lya whose strength varies in time, de- 
pending on the SF history, on the age of the star-forming galaxy , 
and less on its metallicity dSchaererll2003l:lDelgado et aljf20 05). 

The main H and He recombination lines created in the H 11 
re gion surrounding the starburst are also predicted by the models 
of lSchaererl d2003l hereafter S03): for metallicities between 1/50 
Z and solar the strength of Lya varies from EW(Lya)j nt ~ 250- 
360 A at early time after the burst, and declines until zero for 
a burst whereas it reaches an equilibrium value of 60-100 A for 
objects with a constant star formation rate (SFR) after ~ 50-100 
Myr (see also Fig. [TBI). 

In paper II, we showed that the fitting of the observed star- 
forming galaxy cB58 depends only very little on the details of 
the stellar continuum around Lya. Since the objects modeled 
here show stronger Lya emission than cB58, which is dominated 
by absorption, neglecting the detailed shape of the stellar contin- 
uum is even more justified here. Therefore we model the input 
spectrum as a flat continuum plus a Gaussian emission line de- 
scribed by two parameters: 

- the intrinsic equivalent width, EW(Lya)j nt , and 

- the intrinsic full width at half maximum, FWHM(Lya)j nt . 

What "reasonable" values should we adopt for these param- 
eters ? Our first approach, to reduce the number of free param- 
eters in the model, was to test if a unique scenario was con- 
ceivable, i.e. if we could fit all data with the same intrinsic Lya 
spectrum, the differences in the observed spectra would then 
come from radiation transfer effects in the outflowing medium. 
These objects are likely starburst galaxies with a constant star 
formation (SF) his tory as derived from their UV low/medium- 
resolution spectra dNoll et al.ll2004t iMehlert et al.ll2006l) . so we 
fix the intrinsic equivalent width to EW(Lya); nt = 60 - 100 A, 
as derived from the S03 models. We adopted the intrinsic value 
FWHM(Lya)j nt — 100 km s _1 , as it is comparable to the values 
meas ured from the velocity dispersion of Ha and CO lines in 
cB58 dTeplitz et alj|2000b iBaker et alJl2004l). an d the dispersion 
measured in 16 starbursts at z ~ 2 by Erb et al.l d2003l) . 



2.4. Description of the method 

We ran a grid of ~ 500 models with varying physical conditions 
in the shell. The expansion velocity was varied from to 400 
km s _1 in steps of 50 km s _1 , the neutral column density from 
2 x 10 19 to 2 x 10 21 cm" 2 , the dust amount from E(B - V) = to 
0.4 (r fl =0, 0.1, 0.5, 1., 2., 3., and 4.), and the Doppler parameter 
from b = 10 to 200 km s _1 (10, 20, 40, 80, and 200 km s" 1 ). 

In contrast to the four shell parameters, there is no need to 
run a new Monte Carlo simulation each time we want to change 
the input Lya spectrum. This can be made a posteriori with- 
out resorting to any simplifying assumption. To do so we run 
each simulation of the four-dimensional grid with a flat contin- 
uum as input spectrum, e.g. the same number of photons per 
frequency bin, and we memorise this input frequency for each 
photon. Once the simulation is done, we construct output spec- 
tra corresponding to the input frequency bins and we assign a 
different weight to each, in order to reconstruct any input spec- 
trum shape (for example, a flat stellar continuum+a Gaussian 
centered on Lya, or a synthetic starburst spectrum as in paper II. 

From each calculation we derive the integrated spectrum 
(Lya line profile) emerging from the expanding shell. For com- 
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Table 1. Sample of 11 LBG gala xies taken f rom Tapken e t al. (2006) with their observational constraints: ID (col. 1), Lya profile 



type (2), systemic redshift from iNoll etail d2004 except for FDF1267, where z is from T07 (3), the UV (4) and Lya (5) star 
formation rate, the slope of the UV continuum /3 (6), the velocity shift between the LIS lines and Lya, Av(em - abs) (7), the 
observed EW(Lya) Q b s (8) and FWHM(Lya) b. s (9). EWs and FWHM are given in the restframe; we here denote them by "observed" 
for distinction with "intrinsic" or "theoretical" values to be derived later. 



ID 


type 


z 


SFRyv 


SFRLyn- 




Av(em - 


abs) 


EW(Lya) obs 


FWHM(Lycr) obs 








[M G yr'] 


[M 


yr 1 ] 




[km s~ 


'] 


[A] 


[km s~'] 


1267 


C 


2.788 + 0.001 


1.16 ±0.25 


1.49 


±0.08 








129.8 ± 27.41 


235 ± 34 


1337 


A 


3.403 ± 0.004 


27.28 ± 1.15 


2.10 


±0.14 


-2.43 


607 




6.69 ± 0.46 


597 ± 84 


2384 


A 


3.314 ±0.004 


22.74 ± 0.77 


10.8 


±0.27 


-0.55 






83.19 ±3.89 


283 ± 47 


3389 


A 


4.583 ± 0.006 


14.85 ± 2.47 


9.20 


±0.38 








38.82 ± 10.95 


354 ± 70 


4454 


A 


3.085 ± 0.004 


1.98 ±0.49 


2.25 


±0.08 


-2.42 






74.38 ± 11.84 


323 ± 47 


4691 


B 


3.304 ± 0.004 


17.88 ± 0.75 


16.31 


±0.14 


-2.46 






79.44 ± 1.61 


840 ± 115 


5215 


C 


3.148 ±0.004 


26.20 ± 0.80 


9.57 


±0.21 


-1.71 






32.48 ± 1.06 


483 ± 90 


5550 


A 


3.383 ± 0.004 


44.78 ± 1.07 


3.27 


±0.20 


-1.81 


620 




6.36 ± 0.40 


424 ± 85 


5812 


A 


4.995 ± 0.006 


5.24 ± 0.79 


9.60 


±0.18 








153.8 ± 26.6 


226 ± 23 


6557 


A 


4.682 ± 0.006 


13.85 ± 1.39 


3.35 


±0.15 








30.51 ±3.04 


380 ± 135 


7539 


B 


3.287 ± 0.003 


29.87 ± 0.78 


2.45 


±0.46 


-1.74 


80 




6.84 ± 0.46 


1430 ± 230 



parison to the observations, our synthetic spectra are convolved 
with a Gaussian with FWHM=150 km s corresponding to the 
experimental resolution. As shown below, moving in the space of 
these 6 input parameters, we can reproduce the whole diversity 
of observed Lya spectra, ranging from double-peaked profiles to 
broad absorption or asymmetric emission lines. 

Finally, our calculations allow us also to derive the Lya es- 
cape fraction f e . This is computed from 



fe = 



(i) 



where <f>(v) is the intrinsic Lya line profile computed from the 
MC simulation, and f e (v) is the escape fraction in each input 
frequency bin computed from the MC simulation. 

3. Fits of the FORS Deep Field sample 

3.1. Description of the FDF sample 

Our work u s es the FORS Deep Field sample presented by 
Tapk en et alj ([2007, hereafter T07). Therefore, we give a brief 
overview of their work. T07 present the medium-resolution 
spectra (R=2000) of 16 high-redshift galaxies. The target se- 
lection for the objects of T07 was based on the FDF spectro- 
scopic survey (Noll et al. 2004). The FDF spectroscopic survey 
aimed at obtaining low-resolution spectra (R«;20 0) of intrinsi- 
cally bright galaxies with a photometric redshift dBender et al.1 
1200 It iGabasch etaT]|2004l) between z * 1 and 5 with a high 
signal-to-noise ratio (>10). The spectra of 90 galaxies with red- 
shift bet ween 2 and 5 w ere analysed and published electroni- 
cally by dNoll et al.ll2004l) . The deep (up to lOh integration time 
with FORS 1/FORS2) low-resolution spectra allowed them to de- 
rive the redshift with high accuracy and reliability, ad to search 
for even weak signs of AGN activity of the objects. Based on 
the low-resolution spectra, T07 selected starburst galaxies with 
strong Lya emission and/or with bright UV-restframe continuum 
for the follow-up medium-resolution spectroscopy. 

These medium-resolution spectra were obtained with 
FORS 2 at the VLT UT4 using the holographic grisms 1400V 
and 1200R. The spectral resolution of both grisms is R « 2000. 
The spectral range of the 1400V (1200R) grism is about 4500 to 
5800 (5700 to 7300) A. All data were collected in service mode 
using one single MXU mask for each grism. The total integra- 
tion time of the 1400V (1200R) observations is 6.25 h (9.45 h). 



The data were r educed using the MIDAS-based FORS pipeline 
dNoll et al.ll2004T) ). For more details see T07. 

The spectra of all objects of T07 include the Lya- profile. 
However, only eleven Lya profiles have a sufficient SNR (>10), 
which allows a detailed comparison with our theoretical models. 
The properties of this sample are listed in Table Q] Note that 
except stated otherwise all equivalent widths are given in the 
restframe; we denote them by "observed" equivalent widths for 
distinction with "intrinsic" or "theoretical" values to be derived 
later. 

Eight galaxies have redshifts around z»3, while 3 galaxies 
have redshift with 4.5 < z < 5. While only a few Lya profiles of 
our sample show an absorption component (FDF5550), all our 
profiles display an emission component. The equivalent width 
of the emission component range between EW(Lya) D b s = 6 and 
150 A. This Lya equivalent width is measured using the contin- 
uum redwards of the Lya emission line (at x;1300 A). Although 
the majority of LBGs have Lya equivalent widths lower than 
20 A, 8/1 1 of our galaxies have an equivalent width (of the to- 
tal Lya line, including absorption and emission) higher than 20 
A. Therefore 70% of our sample would be detected in a typical 
narrow-band survey, searching for LAEs. 

As described by T07 the profiles show a wide range of 
morphologies. For convenience mostly, we divide the galaxies 
in three groups according to their Lya profile: (A) Lya emit- 
ters with asymmetric profiles: FDF1337, FDF2384, FDF3389, 
FDF4454, FDF5550, FDF5812, and FDF6557, (B) double-peak 
profiles FDF4691 and FDF7539, and (C) asymmetric Lya plus 
a blue bump: FDF5215 and FDF1267. 

3.2. The fitting procedure 

For fits with our synthetic spectra the observed, non-normalised 
spectra were transformed to velocity space using the redshift 
listed in Table Q] If necessary z was adjusted wi thin the error 
bars c ited. We then use the same normalisation as Tapke n^t al.l 
d2007l) to determine EW(Lya)j nt . Finally, we overlay synthetic 
spectra on observed ones and estimate fit qualities. The spectral 
parts we focus on are location of the peak, the shape of the peak 
and the extended wing, knowing that the blue side of the spec- 
trum could be affected by the surrounding IGM. 

The parameters of the best fits, as well as derived param- 
eters from our model like the escape fraction are summarised 
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Table 2. Summary of the best fits derived fro m our model of a spherical expanding shell surrounding a starburst to reproduce a 
sample of 1 1 spectra from Tapke net alj d2006l) . Cols. 1 and 2 are the object ID and Lya profile type respectively. Cols. 3 to 8 give 
the model parameters, col. 9 the derived Lya escape fraction. 



ID 


type 


V cxp [kms '] 


b [km s 1 ] 


Mulcnr 2 ] 




| — ' "1 1 T/T \ TAT 

EW(Lya) mt [A] 


FWHM(LyQ')i n t [km s ] 


fe 


1337 


A 


200 


20 


5 x 10 2U 


1.0 


55 


100 


0.12 


5550 


A 


200 


20 


5 x 10 20 


1.5 


65 


100 


0.05 


2384 


A 


150 


20 


3x 10 19 


1.0 


170 


100 


0.16 


4454 


A 


150 


20 


2x 10 19 


0.5 


100 


150 


0.42 


5812 


A 


150 


20 


2x 10 19 


1.0 


280 


100 


0.16 


3389 


A 


150 


20 


2x 10 19 


0.5 


50 


150 


0.45 


6557 


A 


150 


20 


4x 10 19 


1.0 


70 


100 


0.17 


4691 


B 


10. 


20. 


8x 10 19 


0.0 


80 


1000 


1.0 


7539 


B 


25. 


40. 


5 x 10 20 


0.5 


100 


100 


0.28 


5215 


C 


200 


20 


2 x 10 15 


0.01 


25 


700 


1.0 


5215 


C 


400 


20 


7 x 10 20 


1.0 


120 


100 


0.12 


1267 


C 


50 


20 


2x 10 19 


0.1 


150 


300 


0.64 


1267 


C 


300 


20 


3 x 10 20 


2.0 


500 


100 


0.02 



in Table [2] Multiple entries correspond to multiple solutions of 
similar quality. 

3.3. The asymmetric profiles (group A) 

Among the 1 1 objects, 7 present the characteristic Lye asymmet- 
ric emission line (FDF1337, 5550, 2384, 4454, 5812, 3389 and 
6557). This line shape can be understo od by radiation transfe r 
effects t hrough an ex p andin g medium (IVer hamme et al. 2006). 
Indeed, Tapken et al. (2007) were able to measure a velocity 
shift between the interstellar absorption lines and the Lya emis- 
sion line for two of these objects, FDF1337 and FDF5550, be- 
cause their UV continuum is bright. Both of them present a shift 
of ~ 600 km s _1 (cf. Table [TJ, a clear sign of outflows, which is 
most likely related to 3 times the e xpansion velocity of the shell 
as shown in Verha mme et al.l d2006l) . 

3.3.1. FDF1337 and FDF5550 

To model FDF1337 and FDF5550, we have therefore fixed 
Vexp = 600/3 = 200 km s , and we proceed to adjust the 3 
remaining shell parameters. Our best fits and the corresponding 
parameters are presented in Fig. [TJ Note that both spectra are 
fitted with an intrinsic EW(Lya)j nt ~ 60 A, which corresponds 
to the equilibrium value reached by a galaxy with constant star 
formation after 50-100 Myr (cf. Fig. [TBI. Even if the observed 
EW(Lya) bs is an order of magnitude lower, the intrinsic pre- 
dicted EW(Lyo')i n t seems to be "standard" in this sense. Dust (we 
find T a = 1 .0 for 1337 and 1 .5 for 5550, i.e. E(B-V) ~ 0. 1 -0. 15) 
in a high neutral column density (Nm = 5 x 10 20 cirT 2 ) out- 
flow causes this attenuation. The derived Lya escape fraction is 
f e = 0.052 for FDF5550 and f e = 0.121 forFDF1337. A rather 
small Doppler parameter, b — 20 km s _1 , is derived compared 
to cB58 (£> C B58 = 70km s _1 ) to reproduce the secondary peak on 
the elongated red wing. Note the extension of this red wing over 
~ 1500 km s -1 . 

3.3.2. FDF2384, FDF4454, and FDF5812 

The Lya fits for these objects are shown in Fig. [2] The profiles 
differ from the former by their high EW(Lyo') b s > 70 A, and a 
less extended red wing (~ 600 km s _1 ). This leads to simulated 
neutral column densities an order of magnitude lower and higher 
escape fractions (f e > 0.15, cf. Table |2j. 



To illustrate how well constrained the model parameters are 
we will present in detail the fitting of the object FDF2384. In Fig. 
[3] a fraction of the model grid we use is shown for a fixed value 
b — 20 km s _1 , and for several values of V exp , Nm and t„. The 
central profile of each 3x3 grid illustrates the best fitting profile 
for FDF2384. The observed profile of FDF2384 is overlayed on 
each cell. The overall behaviour of the spectra shown on this Fig. 
can be summarised as follows. 

When V exp increases, multi-peaks appear on the extended 
red wing. Indeed, the location of the second red-peak related to 
"backscattered" photons — photons which are reflected by the 
receding shell through the interior — is at V p ~ 2V exp , so when 
y exp increases, the separation between this peak and the one at 
lower velocity (from the first-red-peak due to scatterings in the 
forthcoming par t of th e shell) increases (for more details, see 
IVerhamme et al.l d2006l) ). 

When Nm increases, the "extension" of the Lya- line in- 
creases — the width between the sharp blue edge and the end 
of the red extended wing. Although very clear, this effect is rel- 
atively modest here, since Nm changes only by a factor 2. This 
effect is easily understood by the increase of the optical depth in 
the medium, forcing Lya to reach frequencies further from line 
center to escape. 

When T fl increases, photons which undergo the highest num- 
ber of scatterings will be destroyed, on average. This enhances 
the peaks made of back-scattered photons, whose escape is eas- 
ier, thanks to this mechanism, than the escape of diffusing pho- 
tons. 

The Doppler parameter is fixed to b — 20km s , because 
higher values broaden the profiles too much, and secondary 
bumps are smoothed. Variations of b are discussed in Sect. |3.6.2l 

The overall shape of the Lya- profile of FDF2384 is smooth, 
with a small bump visible in the extended red wing (at V ~ 400 
km s _I ); its position corresponds to V exp ~ 150 km s _1 . Thanks 
to the "extension" of the line in FDF2384, we can fix Nm ~ (2 — 
4) x 10 19 cirT 2 . Finally, the relative height of the bump compared 
to the main peak will determine correlated values for Nm and 
T a , and two combinations are possible: (N m = 2 x 10 19 and 
T a = 0.5) provides a reasonable fit, but the best fit is obtained 
with (Nm = 3 x 10 19 and t„ = 1.0), as shown in Fig. [2] 

The derived best fit escape fraction is f e ~ 0.16. There are 
big differences of the escape fraction, when t„ increases from 
0.5 to 2, going from f e ~ 40% to f e < 5%; correspondingly the 
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Fig. 1. Lye line profile fits for FDF1337 (left), and FDF550 (right), two of the best constrained cases. The observed spectra are 
shown with the black solid line, model fits as blue dashed curve, the intrinsic (input) profile with the red dotted line. All spectra are 
normalised to unity in the red (positive velocities). The fit parameters are indicated in the Figure and in TableQ] The expanding shell 
model reproduces well the faint and broad asymmetric Lyor emission lines. The secondary peak, or the "bump" in the red extended 
wing is also reproduced for b — 20 km s . Note that the intrinsic EW(Lyar)i nt of these objects is larger than the observed one by 
approximately one order of magnitude. 




Fig. 2. Line profile fits for the three strongest Lya emitters of the sample, presenting all a narrow asymmetric emission line: FDF2384 
(left), FDF4454 (middle), and FDF5812 (right). Same symbols as Fig. inQ] The expansion velocity of the shell is V exp ~ 150 km 
s _1 , similar as for the two precedent objects. The dust content is similar too, but Nui is one order of magnitude lower. The intrinsic 
Lya EW is also larger, particularly for FDF5812 (EW(Lya)mt = 280 A), but these values depend strongly on the continuum 
determination, which is quite uncertain for these faint objects. See text for more details. 



intrinsic strength of the emission line varies from EW; nt ~ 120 
to 170 A. 



We proceed the same way to fit FDF4454 and FDF5812, 
which appear quite similar: the expanding shell has the same 
velocity V exp = 150km s _1 , the column density in the shell is 
an order of magnitude lower than in FDF1337 and FDF5550, 
but the dust amount, t„, remains the same. Since the intrinsic 
EW(Lya)i nt is > EW b s and the latter values are already rel- 
atively large, one obtains quite large intrinsic Lya equivalent 
widths, between ~ 100 and 280 A, for FDF2384, 4454, and 
5812. The highest values require fairly young ages (see Fig. 
[T5t ; however since these objects are quite faint, their continuum 
placement may be uncertain and EWi nt may therefore be overes- 
timated. The Lya escape fraction of FDF4454 is high, f e ~ 0.42, 
because the shell is less dusty. 



3.3.3. FDF3389 and FDF6557 

These two objects are both at z ~ 4.5. They have medium ob- 
served equivalent widths (EW(Lya) b s = 30-40 A), compared 
to the former asymmetric profiles, and rather small extensions 
(~ 600 - 700 km s _1 ). We treat them separately, because their 
spectra seem more complex than a smooth asymmetric emission 
line. They present multi-peaks on the red extension of the line. 
However, they are also more noisy, and these secondary features 
are only twice the noise amplitude. Nevertheless, we assume that 
they are significant, and we derive best-fits, taking these features 
into account (see Fig. |4j. 

The parameters derived from our fitting are close to the 3 
narrow asymmetric cases (V e xp = 150 km s _1 , Nui ~ 2 - 4 x 
10 19 cm' 2 , T fl = 0.5 - 1) except for the intrinsic EW, which 
have "canonical" values — compatible with expectations for 
SFR=const — again, because the observed EW are lower for 




Fig. 3. Grid of predicted Lya line profiles (blue dashed lines) compared to observed spectral line of FDF2384 to illustrate the 
constraints on the fit parameters Nui, Vexp, and r„. All models have been computed with the same Doppler parameter (b = 20 km 
s" 1 ), and the same input spectrum (a flat continuum+a Gaussian emission line with EW; nt =80 Aand FWHMi nt =100 km s" 1 ). Left 
3x3 panel: Variations of V exp (from 100 to 200 km s" 1 from top to bottom line) and r fl (from 0.5 to 2.0 from left to right) for fixed 
Nm = 3 x 10 19 cm" 2 . Right 3x3 panel: Variations of 7V H i ( from 2 x 10 19 cm" 2 to 7V H i = 4 x 10 19 cm" 2 ) and r a (from 0.5 to 2.0 from 
left to right) for fixed V exp = 150 km s" 1 . 



these objects. The escape fraction is high for FDF3389 (f e = 
0.45), and lower for FDF6557 (f e = 0.17) because of a higher 
column density and so a higher dust content. 

3.4. The double-peaked Lya profiles (group B) 

Two of the 1 1 spectra present double-peaked profiles (FDF4691 
and FDF7539). Such line morphologi es are a natura l outcome 
of radiation transfer in a static medium fNeufeldl (1 1990b . since in 
such media the Lya photons can only escape by diffusing into 
the red or blue wings, where the opacity decreases rapidly. 

3.4.1. FDF4691 

This galaxy has a high EW(Lya) obs ~ 80 A. lTapken etaT] d2004l) 
fitted this object wi th a code using a finite element method 
dRichling et alj|200l1) . and proposed 7V H i = 4. X 10 17 cm" 2 and 
b = 60 km s" 1 in an almost static (V exp =12 km s" 1 ) and dust-free 
shell as best-fit parameters. The intrinsic spectrum they use is a 
Gaussian with FWHM=1000 km s" 1 , and no continuum. Using 
the same parameters we can reproduce their fit. 

However, since their code is not suited to high column densi- 
ties we searched for other solutions. We find a fit of better qual- 
ity — the deep gap between the peaks is better reproduced by a 
smaller b, and the red wing is fitted with more accuracy starting 
from an input spectrum with a continuum — with a higher col- 
umn density (Nm = 8 x 10 19 cm" 2 ), and consequently a smaller 
b = 20 km s-fU (see Fig. for the fit and Table [2] for a sum- 
mary of the parameters). Our best fit is obtained with no dust, 



1 Indeed, the location of the peaks were predicted in static media, and 
depend on a combination of Nm and b 



and the fit with r a =0.1 is less good than the one with no dust, 
from which we estimate an upper limit on the dust content of 
T fl 5= 0.1. This is consistent with the fact t hat the Lya and UV 
SFR indicators derived from observations (Tapken et al. 2007) 
are similar. To reproduce the very extended wings of the line, 
the intrinsic Lya emission line has to be very broad. It is char- 
acterised by a very large value of FWHM; nt — 1000 km s" 1 , and 
a "st andard" EW(Lya)i nt = 90 A. Radiation transfer effects are 
inefficient to broaden the line in a medium with such a small col- 
umn density. If interpreted as a result of virial motions, such a 
large FWHM seems, however, unphysical. A hidden A GN may 
be an explanation for the high FWHM, as suggested by iTapkenl 
(120051) . 



A solution to reproduce the observed spectrum with a more 
realistic intrinsic spectrum (FWHMj nt =100km s" 1 andEWi nt =80 
A) is to invoke two contributions from two different media: when 
we sum emergent spectra from two identical shells except for 
the column density (b = 20 km s" 1 , V exp = 10 km s" 1 , no 
dust, EW(Lya) int = 80 A, and nh = 4. x 10 17 cm" 2 for one 
and TVhi = 4. X 10 20 cm" 2 for the other), we are able to re- 
produce a spectrum with narrow peaks close to the center and 
broad wings, starting from a "standard" value for the FWHM, 
FWHM(Lya) int = 100 instead of 1000 km s" 1 (see Fig. ©. 
This could correspond to a physical situation where an initially 
thick shell has been stretched until a hole forms, and the diffuse 
medium in the hole is still opaque enough to imprint radiative 
transfer effects on Lya photons. The surfaces of the thick shell 
and the hole are of equal size in this first model. The parame- 
ters listed for FDF4691 in Table|2]are those of the homogeneous 
single shell model discussed above. 
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velocity (km/s) velocity (km/s) 



Fig. 4. Line profile fits for FDF3389 (left) and FDF6557 (right) showing asymmetric spectra with probable secondary structures on 
the red wing. Same symbols as in Fig.[TJ The parameters derived from our fitting (V e xp = 150 km s" 1 , Nm ~ 2 - 4 x 10 19 cm -2 , 
r a = 0.5 - 1) are similar to the 3 narrow asymmetric cases (cf. Fig. [2]i except for lower intrinsic EW, which here show values 
compatible with expectations for constant SFR over long timescales. 




Fig. 5. Line profile fits for FDF4691 (left) and FDF7539 (right), the two double peaked profiles (type B) with static or almost static 
shells (V exp = 10-25 km s" 1 ). The peak separation and the observed EW(Lyo') bs are different for these objects, and so are the other 
fitting parameters. FDF4691 is the only object for which a very broad input spectrum is derived from the modelling (FWHM = 1000 
km s" 1 instead of ~ 100 for all other objects), maybe a signature of a hidden AGN. 



3.4.2. FDF7539 

For this object the velocity shift between the LIS in absorption 
and Lya- in emission was measured: Av(em - abs) = 80 km s" 1 , 
which implies a shell velocity V exp = 80/3 ~ 27 km s , i.e. 
almost static as in the case of FDF4691, or at maximum ;$ 40 
km s _1 , in case of a low column density. Indeed, the spectrum is 
also double-peaked as for FDF4691. 

The large peaks separation (V p + 500 km s" 1 , larger than 
for FDF4691) implies a high column density. Presumably, the 
rather low observed EW(Lyor) b s also implies the presence of 
dust in the shell. Indeed, the best fit shown in Fig. [5] has a high 
column density (7Y H i = 5 x 10 2l) cm -2 ), and dust (j a = 0.3). 
It is compatible with the canonical value for the intrinsic Lya 
spectrum. The resulting escape fraction is f e = 0.28. 

For comparison, Tap ken et~aT] d2007l) proposed a fit of sim- 
ilar quality for this object, but the velocity of the shell they de- 
rive from their modelling is high (V exp = 190 km s" 1 ), which is 
in contradiction with the observed small velocity shift between 
Lya and LIS (Av(em - abs) = 80 km s~'), and surprising for a 



double-peaked profile. As their investigation is restricted to low 
column densities, the only solution they have to produce sepa- 
rated and broad peaks is with high values of V exp and b, and a 
huge intrinsic FWHM (~ 1900 km s _1 ). 

3.5. Other Lya profiles (group C) 
3.5.1. FDF1267 

This object presents an asymmetric emission peak plus a 
bump on the blue side of this peak with a strong observed 
EW(Lya) b s = 129 A. The Lya profile, shown in Fig. [7j can 
be fitted by different scenarii. 

On the left panel of Fig. |7j the bump is considered as a blue 
peak in an almost static (V e xp = 50 km s~') shell with a small 
column density (Nm — 2 x 10 19 cm" 2 ) and a small amount of 
dust (r fl = 0.1). Values of b higher than b = 20 km s~'lead to 
too separated peaks. The escape fraction is f e - 0.6352. On the 
right panel, we fit the profile with a fast moving (V exp =300 km 
s" 1 ), dense (N m = 3 x 10 2() cm" 2 ) and dusty(r fl = 2.0) shell, 
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Fig. 8. Fitting FDF5215 with two different scenarii : the bump is either considered as a blue peak in an almost static shell (left), or 
as the first red peak of Lya emission (right). 



leading to fe = 0.02, and a very large intrinsic EW(Lyo')i nt 
500 A. Dust is needed in this configuration to reproduce narrow 
and well separated peaks. Indeed, dust destroys more efficiently 
photons which undergo the highest number of scatterings, i.e. all 
but the backscattered photons, which leads "isolates" and "slims 
down" the peaks. Note that this solution requires an adjustment 
of the sou rce redshift to z ~ 2/ 781 instead of z ~ 2.788 + 0.001 
derived bv lTapken et al.l d2007l) . However, this object is the only 
one for which the redshift determination is only based on Lya-, 
so this poses so far no difficulty. 

Our favoured solution is the "quasi-static shell" picture. 
Indeed, the high EWi nt inferred in the second fit seems unlikely. 
Furthermore, the large observed EW b s of 1267 would imply a 
rather low column density as discussed in Sect. 14.2.11 Finally, 
the SFR values derived from uncorrected UV and Lya fluxes 
(SFR(Lya) > S FR(UV)) indicate a low dust content. An accu- 
rate redshift measurement of FDF1267, independent from Lya-, 
should allow to distinguish between these two solutions. 

3.5.2. FDF5215 

FDF5215, shown in Fig. [8]presents the same spectral shape as 
FDF1267: a small bump on the blue side of the asymmetric 
strong emission, but the noise level is much lower, and this small 



bump has to be taken into account. Again, two different scenarii 
can reproduce the spectral shape. The bump is either considered 
as a blue peak in an almost static shell, or as the first red peak 
of Lya emission. The two solutions, differing by more than 5 
orders of magnitude in Nhi, are listed in Table [2] In passing w e 
note that FDF5215 was also modeled by (Tap ken et alj |2007); 
their set of parameters is similar to our solution at low A^hi, ex- 
cept for a higher value of b. 

None of our fits are really satisfying. The solution at Iow-jVhi 
(left panel) reproduces well the observed profile, but the derived 
column density is very low: at least 4 orders of magnitude lower 
than the rest of the sample. As a consequence, the FWHM of 
the intrinsic Lya emission is huge (FWHM=700 km s ) to re- 
produce a broad profile without efficient broadening due to ra- 
diation transfer in an almost transparent medium. Finally, this 
solution does not reproduce the absorption at V = -1200 km s _1 
discussed below. On the other hand, the solution with high A^hi 
and a standard intrinsic FWHM fits less well. Furthermore, the 
redshift derived from this fit is out of the error bars (z = 3.138 
instead of 3.148 + 0.004). 

The black absorption component, found at V ~ -1200 km 
s in this object (see Fig. [8] left), is unaccounted for in the 
fit with low A^hi- How likely it is that this represents a chance 
alignment of an H i absorber? If fitted separately with a Voigt 
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Fig. 6. Fit of FDF4691 with a "two-phase" model, a low density 
shell with Nhi = 4. X 10 17 crrr 2 and a high density shell with 
A^hi = 4.x 10 20 cirT 2 . The other parameters of the shells are iden- 
tical: b = 60 km s _1 , V exp = 12 km s _1 , no dust, EW(Lya)i„ t = 70 
A. This allows for a more reasonable FWHM(Lya); nt = 100 km 
s^instead of 1000. 



profile the absorption is well described by b ~ 70 km s _1 , and 
A^m( abs) ~ 2x 10 17 cm" 2 . Using the column density distribution 
from lMisawa et al] d2007l) we find that the probability to find an 
absorber with > A^Hi(abs) in a velocity interval of say 4000 km 
s is ~ 5 %. It seems thus more likely that this feature is related 
to the galaxy. We conclude that none of our solutions is clearly 
favoured, and further observations are needed to help constrain- 
ing the models for FDF5215. 

3.6. Uncertainties and degeneracies in Lya fits 

3.6.1. Uncertainties on each parameter 

The parameters listed in Table [2] correspond to best-fits deter- 
mined by eye, without resorting to minimisation techniques. We 
now attempt to indicate the approximate uncertainties of the de- 
rived parameters. 

The characteristic uncertainty on the expansion velocity is 
estimated to be ~ 50 km s _1 , which is the step in velocity in our 
grid of models. The sampling in r a , the dust absorption optical 
depth, is not linear (we have assumed r a values of 0, 0.1, 0.5, 
1., 2., 3., and 4. for our grid), but it was refined when necessary. 
From the line fits the characteristic uncertainty on r a and on the 
neutral column density A^hi are estimated as +50%. We refer to 
Fig. [3] to give an idea of the uncertainties on A^hi and T a . Other 
comparisons, e.g. with SED fits and by imposing a consistency 
between different SFR indicators, indicate that the uncertainty 
on the extinction may be somewhat larger, up to a factor ~ 2 in 
some cases (see Sect. 14.8b . 

The Doppler parameter is set to the default value b — 20 
km s _1 , which allows to reproduce the relative narrow peaks ob- 
served , and secondary features on the extended red wing of the 
profiles. Models were computed for other values of b (10, 20, 40, 
80, and 200 km s ), but did not lead to better fits. We estimate 
the uncertainty on this parameter around 50%. 

The uncertainty on the intrinsic equivalent width can be 
fairly large, and it mostly depends on the determination of the 
continuum. Indeed, for strong Lya- emitters, like FDF5812, the 
continuum is so weak that it is poorly constrained, and the un- 
certainty on the continuum level is around 20%. As already men- 



Fig. 9. Dependence of the Lya line profile on b for typical shell 
parameters. Note the shift of the red peak with increasing b, and 
the emergence of a blue counterpart. The escape fraction also 
increases with b. 

tioned, all the spectra presented a bove have been nor malized to 
the same level as determined by Tap ken et alj ([2007 ) to derive 
observed Lya EW b s . However, choosing the continuum level 
by eye we may also obtain acceptable solutions with lower in- 
trinsic EW(Lya); nt ~ 80A for the three stronger Lya emitters 
(FDF2384, FDF4454 and FDF5812), in better agreement with a 
scenario of a constant star formation over -100 Myr. Only for 
FDF1267 do all solutions seem to imply a fairly high intrinsic 
Lya equivalent width. 

The intrinsic FWHM was set to FWHM=100 km s~'as a de- 
fault value, and good fits were obtained for almost all spectra; 
exceptions are the very low A^hi solution for FDF5215 and the 
very broad double-peaked profile of FDF469 1 (except if we con- 
sider a shell with two components as described before). 

3.6.2. Degeneracies 

Although degeneracies affect in principle our profile fits, it turns 
out that asymmetric line profiles provide fairly undegenerate so- 
lutions. We briefly describe how the influence of the main pa- 
rameters (b, Nnh Vexp, FWHM, and t„) can be discerned. 

Doppler parameter b: After several attempts to fit the data 
with large values b (> 50 km s _1 ), we adopted a typical value 
b = 20 km s _1 . The Doppler parameter has a complex influence 
on the Lya profile: small b lead to asymmetric emission lines on 
which the potential multi-peaks due to a high expansion veloc- 
ity would be visible (see Fig. |3). Large values of b (> 50 km 
s _1 ) lead to a smoothed red peak (the multi-peaks are not visi- 
ble any more), whose location is redshifted. Furthermore a blue 
emission component appears whose strength increases with b, 
reproducing a kind of "double-peaked profile" like in static me- 
dia, but with asymmetric peaks — the two sides of each peak 
don't have the same slope (see Fig.|9]l — even in shells with high 
velocities (V e xp > 200 km s _1 ). 

In our sample, the spectra from group A (asymmetric pro- 
files) do not show blue components, and faint multi-peaks (bet- 
ter said bumps on the extended wing) may be visible. Therefore 
small values of b (= 20 km s~') are required to fit our spectra. 
On the contrary, other spectral types (B and C) may be fitted with 
larger values of b. 

Column density A^hi: With increasing neutral hydrogen col- 
umn density, the blue edge of the Lya emission is progressively 
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redshifted with respect to the systemic galaxy redshift. Thus, if 
Z is known accurately enough, the neutral column density in the 
shell is well constrained. The full width of the line increases with 
increasing Nm too. Indeed, it is impossible to fit narrow lines 
(FWHM < 500 km s -1 ) with high column densities (N m > 10 20 
cm -2 ), and extended lines (FWHM > 500 km s~') with low col- 
umn densities (Nm < 10 20 cm -2 ), if the expanding shell model 
applies (cf. Sect. l4.6l and Fig. [13). 

Expansion velocity: The overall shape of the line profile — 
presence of secondary peaks or not — constrains the velocity of 
the shell: fast moving shells (V e xp > 200 km s~') lead to multi- 
peaks in the Lya profile, whereas static or almost static shells 
lead to double-peaked profiles with symmetrical peaks — the 
two sides of each peak have the same slope. 

Dust content: Finally, the dust content is adjusted to fit the 
peak width and the relative height of the bumps - if any - com- 
pared to the main peak. 

In conclusion, few degeneracies appear in the modeling of 
asymmetric Lya line profiles (group A). Thanks to the location 
of the blue edge and to the full width of the line Nm is well 
constrained for asymmetric spectra. Furthermore, the global line 
shape (one single peak) implies low values of b (< 50 km s _1 ) 
and y exp (< 250 km s _1 ). On the contrary, the 4 spectra with more 
comp lex profiles ( 1 267, 5 215, 4691 and 7539) present degenera- 
cies. [Tapken et al. (2007) proposed Lya fits for 3 of them (4691, 
5215, 7539). We can reproduce their fits, but propose fits with 
other sets of paramete rs, as our code allow s for higher column 
densities than the code Tapk en et alJ J2007) used. Automated fit- 
ting methods and a thorough examination of the uncertainties 
and possible degeneracies in Lya line profile fits will be useful 
in the near future, also when larger samples of spectra of suffi- 
cient S/N and resolution become available. 



3.6.3. Possible limitations of the model 

As described and motivated in Sect. [2] our modeling makes some 
simplifying assumptions, including in particular geometry and 
the homogeneity of the shell. How far these assumptions would 
alter our results is presently unclear and remains to be explored 
in the future. 

Inhomogeneous/clum py geometries have e.g. been explored 
by lHansen & Ohl ((2006); the line profiles obtained from such 
models do not seem to change significantly. However, how much 
our model parameters would be modified remains to be exam- 
ined. For the time being it seems clear that few if any cases are 
know, where clumpy geometries w ould favour Lya transmission 
with respect to the continuum (cf. iNeufeldl ll99lHHansen&Ohl 
2006). This can e.g. be co ncluded from the comparison of Ha 
and Lya in local starbursts dAtek et alj 2008). and from the com- 
parison of UV and Lya SFR indicators. Indications for one pos- 
sible case of such a Lya boosting have be en found among 4 
objects analysed by Fin kelstein et alJ d2007al) . Other geometries 
have e.g. been considered in paper II, where deviations from the 
constant velocity shell have been necessary for the analysis of 
cB58. 

The effect of the intergalactic medium (IGM) has been ne- 
glected in our approach. Even if the IGM is almost fully ionised 
at z ~ 3, the redshift of the bulk of our objects, the effect of 
the intervening Lya forest corresponds statistically to a trans- 
mission of ~ 70 and 40 % (r e ff ~ 0.3-1) between z ~ 3 and 4 
(Faucher-Gigu ereet alj|2007). In our modeling we find no need 
to account for such an IGM reduction within ~ 1000-2000 km 
s _1 of the Lya line; not even in the two highest redshift (z ~ 4.7- 
5) objects. No individual Lya forest absorption components are 
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Fig. 10. Comparison of gas extinction, E(B - V), and H i column 
density for the FDF objects (blue triangles), cB58 (blue error 
bar), local starbursts from Calzetti (2007, private communica- 
tion, red cro sses), and measurem ents from the nearby starbursts 
analysed by iKunth et"aT1 (1 19981) (green squares). The multiple 
solutions for the FDF objects are also plotted except for low Nm 
solution of FDF5215; for FDF4691 an arbitrary upper limit of 
E(B — V) < 0.01 is adopted in the plot. The mean Galactic rela- 
tion N H /E(B - V) = 5.8 x 10 21 cirr 2 mag 1 from lBohlin et al.1 
(1 19781) is indicated by the dotted line. 

found in this interval; furthermore for most objects, except pos- 
sibly FDF 1337, 5550, 7539, and maybe also 3389, the predicted 
continuum flux blueward of Lya agrees within the uncertainties 
with the observed continuum. In any case, less importance has 
been given to the line fits on the blue side of Lya. Also, in a 
detailed analysis of t he Lya forest along th e line of sight of the 
z ~ 2.7 LBG cB58, ISavagho et ail (|2002|) found no indication 
for neutral gas within ~ 4000 km s of the systemic velocity of 
the galaxy. From these considerations we conclude that our Lya 
line fits are probably unaffected by additional matter beyond the 
expanding shell included in our models. 

4. Discussion and implications from our model fits 

We now discuss the values of the parameters determined for the 
1 1 objects, possible correlations among them, and we compare 
them with other measurements from the literature. 

4.1. The neutral column density 

The neutral column density we derive from the fitting of 1 1 Lya 
spectra from the FORS Deep Field is radial — along a line of 
sight, from the center of the shell to the end of the simulation vol- 
ume — and ranges over more than an order of magnitude, from 
Nm = 2. x 10 19 to Nm = 5. x 10 20 cm~ 2 . However, the majority of 
objects have a small column density (8/1 1 have Nm < 1 . x 10 20 
cm -2 ). 

How do our Nm determinations from Lya fitting compare 
with other Nm determinations in starbursts? Carrying out such 
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a comparison is difficult for LBGs, since H i column densities 
are usually not measured. Even for nearby starbursts the avail- 
able data is scarce. Calzetti (2007, private communication) has 
kindly determined Nm for us from published 21cm RC3 radio 
observations and assuming sizes given by D25 B F° r compar- 
ison we have also compiled Nm and E(B - V) measurements 
from the s mall sample of near by starbursts observed in the Lya 
region by iKunth et all (Il998l) . These comparison samples are 
plotted in Fig. [10] For S BS 0335-052 we have added a second 
point, ado pting Nm from|Thuan & Izotovl ([1997) and the extinc- 
tion from lAtek et al.l d2008l) . Similarly two points a re shown for 
IRAS08339+65 17 using the extinction compiled bv lKunfh etall 
( 11998b and the one measured by Atek et al. 

Despite some overlap, Nm is lower in our objects than the 
column density observed by Calzetti and collaborators, and more 
similar to the small sample of Kunth et al. An attempt to ex- 
plain this can be the different ways of determining Nm- for the 
Calzetti sample the determination of the neutral column density 
was achieved by radio observations of the whole galaxy, whereas 
in the case of Kunth et al., Nm is derived from Voigt fitting of the 
Lya profile, so it only takes into account the neutral gas which 
influences Lya radiation tra nsfer, even if a Voigt f itting may lead 
to an underestimate of Nm (Verhamme et al. 2006, see Sect. 4). 
This may explain why the determination from Kunth et al. is 
closer to our values than those of Calzetti et al. In any case, to 
compare our column densities with those of Calzetti one needs 
to increase our Nm values typically by a factor ~ 2 to convert 
the radial shell column density to a total one. 

The range of Nm found for the FDF objects is also compat- 
ible with our confirmation of the neutral column density of the 
gravitationally lensed z ~ 3 LBG MS1215-cB58 (cB58, shown 
as the blue cross) that we fitted previously (paper II). Indeed, 
Nm(cB5S) ~ 7. x 10 20 crrT 2 is slightly higher than N m of the 
FDF Lya emitters, as expected for a Lya spectrum in absorp- 
tion. 



4.2. Lya equivalent widths 
4.2.1. Observed EWs 

The observed Lya EWs range from 6 to 150 A in the rest frame. 
Does this range reflect intrinsic differences, or is it somehow 
related to the physical conditions of the ISM in which Lya ra- 
diation transfer takes place ? We examined how EW D b s corre- 
lates with other parameters, but no clear correlation is seen. We 
found a trend in EW D b s with respect to the neutral column den- 
sity in the shell (see Fig. [TTJ: EW obs seems to decrease with Nm, 
at least for the asymmetric profiles (filled circles). In fact, the 
objects with a low EW obs (< 10 A, i.e. FDF1337, 5550, and 
7539) can only be fitted with a high value of Nm, since their 
profiles are very broad. On the other hand, narrow lines with 
large EW b s can only fitted with small values of A^hi- There may 
be three exceptions to this trend, FDF4691, 5215, and 1267. The 
double-peaked profile of 4691 seems peculiar, as it is static and 
dust-free, as suggested by the SFR(UV) and SFR(Lya) which 
are almost identical. We imagine that this trend breaks down in 
dust-free media, or in media with a very small amount of dust, 
because dust is needed to absorb radiation and decrease the in- 
trinsic EW; nt value. FDF1267 and 5215 are peculiar/degenerate 
for the reasons discussed above (Sect. [3]). 



2 The radii obtained in this way reach from ~ 5 to 50 kpc, with a 
median around 17 kpc. 
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Fig. 11. Observed Lya EW b s versus derived radial H 1 column 
density. The three objects with the lowest EW D b s are artificially 
displaced by a small amount for illustration purposes. We may 
see an anticorrelation between the observed Lya EW Q b s and the 
neutral column density in the shell, for the 1 1 objects from the 
FDF, if we eliminate the high-A^i solutions for 1267, and the 
very Iow-A^hi solution for 5215. The filled circles are objects 
presenting asymmetric profiles (type A), the open circles are the 
others (type B, C). The quasi-static object 4691 with a double- 
peaked profile stands out from this trend for unknown reasons. 



An anti-correlation of EW„b s vs Nm can be understood by 
radiation transfer effects if the intrinsic EWj nt is approximately 
constant. When Nm increases, the path length of Lya photons 
increases, and so does their chance to be absorbed by dust: the 
Lya escape fraction decreases with increasing Nm, as mentioned 
above. Furthermore, if we assume a constant dust-to-gas ratio, an 
increase in Nm naturally leads to an increase of the dust quantity 
(the optical depth). These two effects explain the decrease of the 
observed Lya EW with increasing Nm from a theoretical point 
of view. 

4.2.2. Intrinsic EWs 

Three objects (2384, 5812, 1267) have clearly very large intrin- 
sic EW(Lya)j nt > 100 A, for one (5215) the two solutions give 
quite disparate results, and the remaining 7 objects have all in- 
trinsic EWim of ~ 50-100 A. Taking the uncertainties in the con- 
tinuum placement into account (cf. above) we consider that this 
latter group (7 of 1 1 objects) have intrinsic EW(Lya)j nt compat- 
ible with expectations for star-forming galaxies with a constant 
star formation history over periods 10-100 Myr, as seen in 
Fig-El The strength of Lya in three high EW Q b s objects requires 
younger ages, irrespectively of their star formation history. 

4.3. Dust extinction 

For the first time we have derived here constraints on the dust 
content of galaxies using the Lya profile only. It is therefore of 



Verhamme, Schaerer, Atek, Tapken: 3D Lya radiation transfer. III. Gas and stellar properties of z ~ 3 LBGs 



13 




0.1 0.2 0.3 

E(B-V) - (Hyperz or /S) [mag] 



Fig. 12. Comparison of the extinction E(B - V) determined from 
the Lya profile fits versus other methods for objects with suffi- 
cient photometry and/or measured B slopes. Multiple solutions 
from Lya fits are included. Hyperz SED fits using constant SFR 
models (red filled triangles) or arbitrary SF histories (blue open 
circles). Green squares show values of E(B - V) gas determined 
from the /3-slope. The dotted line shows the one-to-one relation, 
the solid line t he relation E(B - V) * = 0A4E(B - V)s as found 
empirically (cf. lCalzetti et al .1120001) . 

interest to examine how this determination compares with other 
methods. 

Using a version of the SED fit and photo metric redshift code 
Hyperz described in lSchaerer & Pellol 00051) and the UBgRIJKs 
photometry published by iHeidt et all d2003l) . we have modeled 
the SED of our obj ects, assuming the Calzetti attenuation law 
dCalzetti et al.ll2000l) . Three objects, 5812, 3389, and 6557, have 
insufficient data (photometry in 3 or less bands) which does not 
allow meaningful SED fits. 2384 is also excluded, since it ap- 
pears to be a multiple source, where the Lya emission is clearly 
displaced from the continuum. Results for the remaining objects 
are shown in Fig. Q~2] where we compare the extinction derived 
assuming models with constant star formation (red filled trian- 
gles) or exponentially decreasing SF histories (blue open circles) 
described by the Bruzual & Chariot templates with E(B - V) de- 
rived from our Lya line fits. Note that these values, denoted here 
as E(B - V)*(Hyperz), measure the extinction suffered by the 
stars, whereas E(B - V)(Lyo') measures that of the Lya emitting 
gas. The two may differ, as e.g. known to hold empirically for 
local starbursts between the stella r extinction and the one mea- 
sured from the Balmer decrement dCalzetti et al. 2000). 

Figure Q~2] shows a good correlation between the different 
extinction measures, especially when different SF histories are 
allowed for. Indeed, for the bulk of the objects the extinc- 
tion derived from Lye profile fitting is between the gas extinc- 
tion expected from Calzetti's empirical relation, E(B - V) gas = 
1 /0ME(B - V)*, and a somewhat lower value of E(B - V) gas . 

We can also estimate the extinction from the UV slope. 
Excluding again the multiple source 2384, the observed UV 
slopes B show basically two groups, whose extinction, estimated 



following lCalzetti et al.ld2000l) . is E(B-V)* ~ 0.-0.02 and 0.16- 
0.18 shown by the green squares in Fig. [12] These values cover 
a similar range as E(B - V)*(Hyperz), although with a poor cor- 
relation. This could be due to the statistical nature of the under- 
lying correlation between attenuation and B. We conclude that 
overall Lya line profile fits allow us to obtain quite consistent 
extinction values compared to broad band photometry fits of the 
individual objects. Our derived extinction values, corresponding 
to E(B - V) ~ 0-0.2 ar e also in good agreem ent with the val- 
ues found for LBGs by IShaplev et al.l d2003l) and others (e.g. 
iPapovich et ai]|200ll) . 

4.4. Gas to dust ratio 

How do the gas-to-dust ratios obtained from our line profile fits 
compare with other values observed in starbursts? To address 
this we turn again to our comparison samples shown in Fig.fTUl 

Compared to the Galactic average of Nhi /E(B - V) - 5.8 x 
10 21 ctrr 2 mag" 1 and its ~ 30% scatter (cf. Bo hlin et al .11 1978b , 
our results show somewhat lower gas-to-dust ratios and appar- 
ently a larger scatter. The scatter is, however, similar to the 
one found among the local starbursts also shown in this Fig. 
Petti ni et"aT] (|2000|) have already noted this difference in gas- 
to-dust ratio for cB58, which they suggest could indicate that a 
significant fraction of the gas is not in atomic form, i.e. is either 
ionised and/or in molecular hydrogen. 

The Calzetti objects have a median of Nm/E(B - V) = 
7.9 X 10 21 similar to the "classical" Galactic value but show- 
ing a wide dispersion. Approximately half of our objects show 
Nhi/E{B - V) values lower than those of Calzetti's local star- 
bursts. For the other half of the sample, and for the LBG cB58 
modeled in paper II, the values of the gas-to-dust ratio overlap 
with those for the Calzetti sample. In any case the comparison 
may be hampered for the reasons affecting also the Nm compar- 
ison (cf. Sect. 14. U . Overall we conclude that the gas-to-dust ra- 
tios determined purely from Lya line profile fitting yield values 
reasonably cons istent with local st arbursts or somewhat lower. 
For comparison I Vladilo et al.l (120071) find higher gas-to-dust ra- 
tios in DLAs compared to the Galactic value. 

If some LBGs show truly lower gas-to-dust ratios, this may 
be due to a higher degree of metal enrichment in the outflow- 
ing H i gas, and/or an overall smaller fraction of neutral gas be- 
ing "polluted" in nearby starbursts, and/or a smaller fraction of 
atomic hydrogen. However, the column densities found in LBGs 
are also consistently lower than in local starbursts (cf. Fig. [TOl. 
which may indicate a different "regime". Other, independent de- 
terminations of the gas-to-dust ratio in LBGs and other starbursts 
and more detailed examinations would be required to clarify 
these issues and to understand the possible physical causes for 
these apparent differences. 

4.5. Velocity of the outflow 

Assuming a galaxy-scale outflow surrounding our 11 objects 
(the relevance of this model is discussed in Sect. |2|, we have 
derived the expansion velocity V exp of the expanding shell from 
the Lya spectral shape for 8 objects, and deduced it for the 3 
other objects (1337, 5550, 7539) from the measurement of a 
shift between the LIS absorption lines and the Lya emission. 
Overall 9 objects out of 1 1 have shell velocities around 150-200 
km s" 1 , and 2 objects present almost static shells. No correlation 
between the expansion velocity and other parameters is found. 
We now briefly discuss the high velocity outflows and the few 
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nearly static cases. Beforehand we can already mention that we 
have no explanation for the causes leading to low ISM veloc- 
ities or allowing to distinguishing these objects from the more 
common cases showing outflows. 

4.5.1. High velocity objects 

All objects presenting an asymmetric emission line are repro- 
duced with V exp ~ 150 - 200 k m s" 1 . This velocity ra nge is very 
similar to the determination of Shapl evet al.l (|2003), from the 
blueshift of LIS absorption lines with respect to stellar lines in 
a sample of ~ 800 LBGs at redshift z ~ 3, as well as in cB58, 
where the outflow velocity is estimated to V exp = 255 km s" 1 
dPettini et alj|2002l) . 

Two objects with a more peculiar spectral shape, FDF 1267 
and 5215, can also be fitted with higher expansion velocities 
(Vexp ~ 400 km s" 1 ), a high Nm, and high dust content. However, 
1267 has such a large observed equivalent width (EW b s = 129 
A) that it may rather have a small column density to fit better in 
the plot showing a correlation between the observed Lya EW b s 
and Nhi (see Fig.fTTI). In passing we note that this peculiar spec- 
tral shape seems als o to be found i n ob servations of an LBG at 
redshift z ~ 3.7 bv lVanzella et alj d2008l) in the GOODS-South 
field. 



4.5.2. Low velocity objects 

The two double-peaked spectra (4691 and 7539) are charac- 
terised by a static (or almost static, V exp < 25 km s" 1 ) surround- 
ing shell. In the case of 7539, the velocity shift between the LIS 
absorption lines and Lya is even measured (Av(em - abs) = 80 
km s" 1 ), so the shell velocity is here an observational constraint 
(V eX p ~ Av(em - abs)/3). This peculiar spectral shape (dou- 
ble peaks) is predicted by theory, arising from L ya resonant 
scattering through static H i media (Neufeld 1990), but surpris- 
ingly, observed double-peaked profiles are in general not inter- 
pret ed as a signature of Lya radiation transfer through static me- 
dia dFosburv et alj|2003l Ichristensen et al.l|2004t IWilman et all 
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2005: IVenemans et alj|2005t IVanzella et al.ll2008l) . even if they 
appear much less common than the asymmetric emission ob- 
served in all high-z LAEs. In the static case the separation of the 
peaks AA is not only related to the thermal velocity of the H i 
gas, but depends also on the H i column towards the Lya source, 
A^i- For a homogeneous slab one has: 



AA R 



2A (-x„b/c) 



2.49 x 10" 



12.85 km s" 



W (cm- 2 ) 



A 



(2) 



where c is the light speed, Aq is the restframe wavelength of 

Lya, b - JV^ + Vf urb is the Doppler parameter, and x p = 

0.88(q to)^ 3 is the location of the peaks in u nits of the Doppler 
width dNeufeldlll99llDiikstra et al.ll2006bl) . For example, for 
b = 12.85 km s~'and N m = 6.4 x 10 19 cm" 2 one obtains 
AA KSt ~ 1 A. If applied to t he z ~ 3.65 double peaked Lya 
object of Vanz ella et al.l d2008l) . the observed velocity shift of 1 3 
A would indicate Nm ~ 1.4 X 10 21 cm" 2 for b = 12.85 km s" 1 . 
Establishing accurate enough galaxy redshifts for these objects 
is important to be able to assert if one is truly dealing with a 
nearly static case (in which case zero velocity is between the 
two Lya peaks), outflows (with both Lya peaks redshifted), or 
other situations. 
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Fig. 13. Observed Lya FWHM versus neutral column den- 
sity, showing a tentative correlation between the observed Lya 
FWHM and the H i column density in the shell, for the FDF ob- 
jects. The filled circles are objects with asymmetric Lya profiles, 
the open circles are the others. The objects with double-peaked 
profiles, FDF 4691 and FDF 7539, are clearly distinct, showing 
the highest FWHM. 



4.6. Observed Lya FWHM 

A possible correlation may be found between the observed 
FWHM and the neutral column density of the expanding shell, 
as shown in Fig. Q~3] If real, it may be used to estimate the 
H i column density in starbursts from a simple measurement of 
FWHM(Lya). Such a correlation can easily be explained by ra- 
diation transfer effects: the Lya optical depth increases with Nm 
so that Lya photons have to diffuse further in the wings to be 
able to escape the medium, which naturally broadens the Lya 
red wing. 

Another point worth noticing about FWHM is that we were 
able to fit all observed Lya profiles with an intrinsic FWHM of 
~ 100 km s" 1 , except for 4691, for which we have to start with 
an already very extended intrinsic Lya line to reproduce the very 
broad double-peaked profile. Otherwise, more exotic scenarii, 
like a shell with a hole, have to be invoked. 

Finally, the two correlations presented above (FWHM D b s 
vs Nm and EW b s vs Nm, Figs. QTJ and [TBI explain the ob- 
served anti-correlation shown by T07 between FWHM b s and 
EW bs- Indeed, they are both related to the neutral column den- 
sity A^hi which surrounds the starburst. The EW D b s increases 
with decreasing Nm and the FWHM D b s increases with increas- 
ing A^hi- Therefore, we should not observe objects with a large 
EW b s and very broad lines, what is also confirmed by the ob- 
served FWHM b s of LAEs which are always below 500 km s" 1 
dRhoads et al.ll2003l:lDawson et al.ll2004tlVenemans et al.ll2 004). 
If the tentative correlation between FWHM b s and Nm really 
holds it would imply a maximum column density of Nm ~ 
(2 -4)x 10 20 cm" 2 in LAEs. 
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Fig. 14. Lya escape fraction versus dust extinction in the gas 
for the LBG (circles) and local starbursts (crosses). We find a 
clear correlation between the Lya escape fraction and the dust 
amount in the shell, for the 1 1 objects from the FDF. The filled 
circles stand for objects with asymmetric profiles, open circles 
the remaining ones. The solid line represents the continuum at- 
tenuation, fe cont = exp(-T n ) » exp(-10 x E(B - V)), the dashed 
line the fit proposed in Eq. [3] The crosses and the upper limit 
are the inte grated escape fract ions from a sample of 6 local star- 
bursts from lAtek et al.1 12008) plotted as a function of E(B - V) 
measured from the Balmer decrement. 

4. 7. Escape fraction 

The Lya line escape fractions derived for our FDF objects range 
from 100% for the dust-free object FDF4691 to ~ 2% for ob- 
jects with the highest extinction. From our modeling we find 
that the main parameter determining the escape fraction is the 
dust amount in the shell. Furthermore, from all our modeling re- 
sults already discussed above, it is quite clear that no single value 
of the Lya escape fraction is expected from LBGs and LAEs, in 
contrast to simplifyi ng assumptions made in some models (e.g. 
iLe Delliou etal]|2006b . 

Our model grids predict the following behaviour for the es- 
cape fraction: 

- f e increases with increasing V exp , because the Lya optical 
depth decreases in a fast moving medium compared to the 
static case. Therefore the mean path of Lya photons in the 
medium decreases, so their chance to be absorbed by dust 
too. 

- f e decreases with increasing Nhi, because the Lya optical 
depth is proportional to Nm- 

- f e decreases with increasing E(B - V), obviously since with 
a larger number of absorbers in the medium, Lya photons 
increase their chance to interact with them. 

However, the only clear correlation we found in the data between 
f e and other parameters is with E(B - V) (see Fig. [T4b. No cor- 
relation is found in particular with V exp (our objects cover prob- 
ably a too small range in velocity range), and with Nhi (the vari- 
ation in the dust-to-gas ratio is probably more dominant). On 



Fig. [14] we see also that, as expected, Lya photons are more 
attenuated by dust than the continuum; indeed all LBG data 
points (circles) are located below the solid line corresponding to 
/ e (cont) = exp(-T fl ). The reason is that multiple resonant scat- 
tering of Lya photons increases their path through the medium, 
and hence their chance to be absorbed by dust, compare to the 
continuum. Note that other dependences of f e on E(B - V) could 
be expected with other geometries. For example in clumpy me- 
dia, the reflection of Lya on the clump surfaces could ease the 
Lya transmission fo r the same amount of dust (Neufeld 199U 
lHansen & Ohll2006l) . 

We propose a fit to predict the escape fraction of Lya photons 
knowing the dust extinction (dashed curve on Fig. FBI: 

f e = i - 7 - 7I x^>. (3) 

Note that E(B-V) in the precedent formulae is the extinction in 
the gas, which may be different from the extinction of the stars 
(ICalzetti et al. 2000), as already mentioned above. Interestingly, 
the two static objects are also fitted by this formulae, which illus- 
trates that dust is really the dominant parameter which governs 
the Lya escape in our objects. One of these (4691) is dust-free, 
so its escape fraction is ~ 1, but in the other object (7539), 30% 
of the Lya photons escape the medium. For the same extinction, 
moving media present an escape fraction of 40-45%, which is 
coherent with the theoretical prediction that f e increases with 

Vexp- 

Empirical Lya es cape fractions have recently been measured 
bv lAteketall J2008) from imaging for a sample of 6 local star- 
bursts. Their values are compared to our data for LBGs in Fig. 
[T4l For E(B — V)< 0.2, our results are in good agreement with 
three local objects. SBS 0335-052 with an integrated extinction 
of E(B - V) gas ss 0.21 shows no Lya emission, it is a net ab- 
sorber. For larger E(B - V) values, the Lya escape fraction of 
two local starbursts (Haro 1 1 and NGC 6090) are higher than f e 
predicted by our fit to the LBGs studied here. Deviations from a 
simple homogeneous shell geometry are the most likely expla- 
nation for this difference. This will be testable through detailed 
modeling both of the spatially resolved and integrated properties 
of the local objects. 

Since the Lya line flux is more strongly reduced (due to mul- 
tiple scattering effects) than the adjacent continuum, the Lya 
equivalent width depends on the extinction. This phenomenon is 
added to the one already known to result from the extinction dif- 
ference between the gas and the stellar continuum. In principle, 
a measurement of EW(Lya) b s could thus be used to determine 
the extinction, provided the intrinsic equivalent width EWj nt is 
known. Concretely, the fit-relation between f e and E(B - V) gas 
proposed above (Eq. [3) translates to the following behaviour of 
the Lya equivalent width with extinction: 

/ EW rest \ 

\ EWi m / 

* -5.6 E(B - V) gas (4) 

where ^ * 12 an d r =E(B-V)*/E(B-V) gas = 0.44 ac- 
cording to Calzetti ( 200 1 ) . Adopting reasonable values for E W; nt 
(e.g. from Fig. [15), this formula may be used to obtain a crude 
estimate of the extinction in LAEs based on a pure equivalent 
width measurement. An extinction corrected SFR(Lya) value 
can then be obtained from the Lya luminosity using an appropri- 
ate SFR calibration from Fig. [15] consistently with the assumed 
value of EWi nt . 
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Z = 0.02, 0.004, 0.0004 (S03 models) 
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Fig. 15. Temporal evolution of Lya and UV SFR predictions 
from the synthesis models of S03 for three metallicities (Z = 
0.02 = Z in black, 0.004 in red, and 0.0004 in blue) computed 
for instantaneous bursts and/or constant SF. Top: Lya- line lumi- 
nosity in erg s emitted per unit SF rate, assuming a Salpeter 
IMF from 0.1 to 1 00 M Q . The dott ed line shows the "canoni- 
cal" value based on iKennicutlld 19981) and a standard Lya/Ha ra- 
tio. Middle: logarithm of the UV to Lya SFR ratio. The shaded 
area shows allowed range allowed for constant SF models with 
metallicities between 1/50 Z and solar. Bottom: Lya equivalent 
width. 



Fig. 16. Comparison of the SFR values determined from the UV 
continuum and from Lya for our objects. Blue circles indicated 
the "observed", uncorrected SFR values from Table Q] Green 
symbols show the "true" values corrected for dust and transfer 
effects according to Eqs. [5] and [6] The shaded area shows the 
range allowed by synthesis models for the combination of the 
"true" SFR(Lya) and SFR(UV) values taking age effects of con- 
stant SF models into account and allowing for metallicities be- 
tween 1/50 Z Q and solar (cf. Fig.fTBIl. The dashed line indicates 
the one-to-one relation. 



4.8. SFR indicators 

Given our quantitative analysis of Lya- radiation transfer, the de- 
termination of the Lya escape fraction and of the extinction, we 
are now able to examine to what extent Lya and the UV contin- 
uum provide consistent measures of the star formation rate. The 
main results of this exercise are shown in Fig. [16] 

First we note that three of our objects (1267, 4454, 5812) 
show observed, i.e. uncorrected SFR values corresponding to 
SFR(Lya) >SFR(f/V). Such a result need not be inconsis- 
tent; this behaviour is indeed expected for young bursts or 
objects where constant star formation has not yet proceeded 
over long enough timescales, i.e. for timescales ^ 10-100 Myr 
as shown in Fig. Q3] In this case values up to SFR(Lya) ~ 
4xSFR(£/V) can be obtained; this allowed range of values for 
SFR(Lya)/SFR(t/V) for constant SF is shown by the shaded 
region on Fig. [16] Indeed these objects also show among the 
largest EW(Lya) t, s , as expected from Fig. Q3] for relatively 
young, but on average constantly star-forming, objects. Second, 
4691 shows SFR(Lya) * SFR(UV) (observed), indicative of lit- 
tle or no dust, and confirmed by our modeling. Finally the re- 
maining 7 objects show SFR(UV) > SFR(Lya) with UV star 
formation rates up to ~ 14 times larger tha n Lya, a result found 
for the majority of LBGs and LAEs (e.g. Ya mada et al.l 12005; 
lGronwalletailf2007h . 

We now correct these SFR indicators for the effects of dust 
and radiation transfer. One has 

SFR(Lya) true = SFR(Lya)// e , (5) 



with the Lya escape fraction f e . The UV SFR is corrected assum- 
ing the Calzetti law and the extinction E(B - V) derived from the 
Lya line fit, i.e. 

SFR(UV) true = SFR(UV) x io+o^(.B-vrk vv ^ (6) 

with E(B-V)* = 0A4E(B-V). For our objects k vw = k l60l) a 10 
is appropriate. Using E(B - V) derived from our Lya fits, the re- 
sulting "true" SFR values for all our objects are plotted in Fig. 
[I6](green triangles). For the majority of the objects we find that 
their corrected SFR values show much less dispersion between 
Lya and UV based measurements than observed, uncorrected 
values. Three objects (1337, 5215, 7539) seem to require even 
lower Lya escape fractions (i.e. larger extinction) to reduce the 
differences between their SFR indicators further. Adopting Eq. 
[3]and imposing SFR(Lya)=SFR(UV) would imply an extinction 
E(B - V) of (0. 19, 0.07, 0.18) compared to our estimate of (0. 1, 
0.1, 0.05) mag from line fitting. Similarly the extinction seems 
to be overestimated in FDF 5812 and for the high Nhi solution 
of 1267. The maximum extinction allowed to obtain SFR consis- 
tencies for these objects within the shaded region is E(B - V) = 
0.09 and 0.06 for 1267 and 5812 respectively. Most of these "ad- 
justments" are within a factor of 2 uncertainty. Overall we there- 
fore conclude that the results based on our Lya line fits using 
radiation transfer models improve the consistency between UV 
continuum and Lya line based SFR indicators. 
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4.9. Other correlations 

Several correlations have been found earlier between Lya and 
other properties of LBGs from analysis of large spectroscopic 
samples. When groupin g their ~ 1000 LBG spectra according to 
Lya equivalent width, IShaplev et al.l (120031) found the following 
main correlations with decreasing EW(Lya) b S (i.e. from emis- 
sion to absorption): 1) The extinction increases, 2) the strength 
of the low-ionisation interstellar lines increases, 3) the veloc- 
ity shift between these interstellar lines and the Lya peak in- 
creases Av(em - abs), and 4) the dust-corrected SFR increases. 
Corr elations 1, 2, and 4 have also been fo und in other samples 
(e.g. iNoll et all 120041; rtapken et al.ll2007l) . We now examine if 
these correlations also hold for our small sample, and how our 
model may or may not explain them. 

4.9.1 . E(B - V) versus Lya EW 

A decrease of EW(Lya) b s with increasing extinction is natu- 
rally expected from our radiation transfer models, e.g. when the 
intrinsic equivalent width is constant (see Eq. 2]). This seems 
the most simple/natural explanation for the trend observed be- 
tween the average LBG spectra, also given the relative long SF 
timescales observed generally in LBGs, which should imply a 
fairly constant intrinsic EW(Lya)j nt . 

The 1 1 individual LBGs analysed here show no clear trend 
between EW(Lya) b s and E(B - V). Given the relatively large 
scatter in the derived gas-to-dust ratio (see Fig.fTOli it appears that 
the clearest observational trend found is between EW Q b s and Nm 
instead (Fig. ITTb . Of course, if on average the gas-to-dust ratio 
is constant, any correlation between EW Q b s and Nm would auto- 
matically impl y a correlation betwee n E(B-V) and EW(Lya) b s , 
as observed bv lShaolev etai](l2003l) . 

4.9.2. EW(LIS) versus EW(Lya) 

Our model does not make direct predictions for the strengths 
of interstellar absorption lines. For this reason and since most 
of the observed low ionisation interstellar lines are saturated 
(Sha plev et alj|2003l) we are not able to examine quantitatively 

thi s correlation. 

iFerrara & Ricottil (|2006) propose to explain the behaviour 
of EW(LIS) by cold debris whose covering factor decreases 
as a function of time in a dynamical outflow model. EW(LIS) 
is then mostly related with the wind velocity (constrained by 
Av(em - abs)) which may vary with EW(Lya) b s according to 
these authors. Verifying observationally whether and how wind 
velocity truly varies and quantifying the covering factor of the 
cold outflowing gas is therefore important to test this scenario. 

4.9.3. Av(em - abs) versus EW(Lya) 

Shap lev et all (120031) measure velocity shifts between 475 ± 25 
and 795 + 3 km s for their four groups (quartiles) of LBGs 
with decreasing EW(Lya) b s - To verify this trend we have re- 
measured Av(em - abs) and oth er quantities from the composite 
spectra of IShaplev et al.l d2003l) : overall our measurements are 
within lcr, except for group 1 (with the strongest Lya absorp- 
tion) where we obtain Av(em - abs) = 703 km s _1 . Clearly the 
uncertainty cited for this group must be underestimated. Does 
this observed trend imply a systematic variation of the average 
outflow velocity between the spectral groups? For example, if 
true, this could imply an increase of the average wind velocity 
from V exp » 1/3 X Av(em - abs) ~ 160 to 235 km s _1 (or to 
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Fig. 17. Dependence of the Lya line profile on Nm for typical 
shell parameters. Note the clear shift of the peak of the emission 
profile from low velocity (~ 100 km s _1 « 0.5 x V ex p) to twice 
the shell velocity with increasing Nm, typical for Nm ~ 10 20 
cirT 2 . 



265 km s 1 using Shapley's value for group 1) if spherical shell 
models apply as argued earlier. 

As already mentioned above, the expansion velocities de- 
rived for our objects and those measured for 2 of them as well 
as the t ypical reddening val ues, are in good agreement with the 
data of Shapl ev et all (120031) . However, within our small sample 
we do not find any systematic variation of V exp , and clearly the 
reddening variations we find are large compared to the differ- 
ence of the mean extinction between groups 3 and 4 of Shapley 
(where our objects would lie) and compared to the dispersion of 
E(B - V) within these groups. 

Shapley stresses that "none of the correlations with Av(em - 
abs) is as significant as the trends among EW(Lya) b s , EW(LIS), 
and E(B-Vy. For example, when grou ped according to Av(em- 
abs), the objects of Shap lev et al.l (120031) show several trends op- 
posite to those trend found between the groups constructed ac- 
cording to EW(Lya) bJl The reality and significance of this be- 
haviour may thus be questionable. In fact we find two effects, a 
radiation transfer and an "instrumental" one, which could artifi- 
cially lead to such a behaviour. 

First, while for expanding shells with column densities 
A^hi ^ 10 20 cirT 2 the main peak of the Lya profiles is indeed 
redshifted by twice th e expansion velocity (as pointed out by 
IVerhamme et al.l 120061) . the peak emerges at lower velocities 
(~ 1 x Vexp) for lower column densities. This transition is clearly 
illustrated in Fig. [17] If Shapley's spectral groups correspond on 
average to a sequence with increasing Nm (from Lya in emission 
to absorption) and if this transition happens somewhere within 
this sample it would mean that the real spread of outflow veloc- 



3 E.g. Instead of decreasing, E(B - V) increases with increasing 
EW(Lyo') b s between LBG s grouped according to Av(em-abs) (cf. Figs. 
3 and 16 of Shapl ev et alj d2003h ). 
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LBG and LAE samples at z ~3 — 5 
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Fig. 18. Effect of spectral resolution on the Lya line shift mea- 
surement. The dotte d line is the observed spectrum of cB58 from 
iPettini et all d2002l) with a resolution of 58 km s _1 . Red solid 
lines show the same spectrum convolved to resolutions of 100, 
200, 400, and 700 k m s~'. The blue das hed line is the spectrum 
of the first group of Shaplev et al. (2003) normalised to the one 
of cB58 at +1700 km s _I . Note the progressive redshift of the 
Lya peak (observed at +355 km s _1 in the original cB58 spec- 
trum) by up to ~ 200 km s _1 for a resolution of 400 km s _1 , 
typical of the composite spectra of Shaplev ^t al] d2003l) . 



ities would considerably be reduced to V exp from ~ 238 to 265 
km s _1 (or 235 with our measurement for group 1). 

A second, unavoidable effect, also alters the observed veloc- 
ity shift between Lya and the interstellar lines. Indeed, given 
the increasing strength of broad Lya absorption in groups 2 
and 1 compared to the other groups, the relatively low spec- 
tral resolution of the composite spectra affects the measure- 
ment of the Lya line fit introducing a systematic shift. For il- 
lustration we show in Fig. [18] the high-resolution spectrum of 
cB58 - an LBG with Lya- properties characteristic of group 1 
- and the effect lower spectral resolution has on the measure- 
ment ^f_fhe_I^y_a_geak 1 Jnterestingly but not really surprisingly 
(cf. iMas-Hesse et al] [20031) . a limited spectral resolution, such 
as the one of ~ 400 -700 km s _1 of the composite spectra of 
Shapl ev et al.1 (HH, can lead to an overestimate of the true 
Lya line redshift in objects with a steeply rising "underlying" 
spectrum. In other words, due to the appearance of the broad 
Lya absorption wing, the measurement of Av(em - abs) may be 
affected by a systematic shift which increases with decreasing 
EW(Lya) bs i.e. towards the group 1 with the strongest Lya ab- 
sorption. As the exercise in Fig.[18]shows, this shift can be of the 
same order as the apparent velocity differences between different 
LBG subgroups. For the two reasons just discussed we conclude 
that the observed variatio ns between Lya and the IS lines found 
between the subgroups of Shaplev et al. (2003) and possible cor- 
relations with other quantities cannot reflect a simple behaviour 
of the real outflow velocity. More work is needed to clarify the 
possible link between outflow velocities and other properties of 
LBGs. 



Fig. 19. Observed restframe Lya equivalent width of LBGs and 
LAEs as a function of their UV star formation rate (uncor- 
rected for reddening ). The following data is shown: LBGs from 
iTapken et al] d2007) (blue triangles), cB58 IPettini et al] d2002l) 
(blue circle), and the individual objects from I Shaplev et al] 
(200 3]) (magenta cr osses). LAEs o r LAE-cand i dates from 
lYamada etaLl d2005l) (red triangles) lOuchi et al] d2007l) (red 
squares) and from iGronwall et all (1200717 (black circles). The 
data clearly shows the absence of high equivalent width objects 
at high SFR, which we explain as being due to dust/radiation 
transfer effects. For further details see text. 



In any case our radiation transfer calculations show that the 
variations of V exp found between the objects fitted here have a 
relatively minor influence on properties such as EWiLya), the 
Lya escape fraction etc. It appears that as long as global out- 
flows with velocities of several 100 km s _I are present, other 
parameters are more dominant. 



4.9.4. SFR(UV) versus £W(Lya) 

Trends or correlations between EW(Lya) b s and the star for- 
mation rate have been noted by many groups and are found 
not only for LBGs, bu t also for LAEs (e .g. IShaplev et al] 
120031: lAndo et al] 12004 ITapken et all 120071) . A compilation 
of available data is shown in Fig. [19] whe r e we include 
the following data : LBGs from ITapken et al] (120071) . cB58 
IPettini et al j ( 20021) (blue symbols), and the individual objects 
from IShaplev et all d2003l) (m a genta c rosses); LAEs or LAE- 
candidates from I Yamada et al] d2005l): Ouchi et all d2007l) (red 
symbols) and from Gronwall et al] (2007) (black circles). Note 
that Tapken's LBGs include four z ~ 4.5-5 objects, Yamada's 
sample of 198 LAEs spans z ~ 3.3 to 4.8, Ouchi's 84 spectro- 
scopically confirmed LAEs have z ~ 3.1, 3.7 and 5.7. All other 
objects have redshifts have z ~ 2.7-3.2. Note that none of the 
data was corrected for reddening. SFR(UV) is therefore directly 
proportional to the observed UV magnitude of the objects. 

The observations summarised in Fig.[T9lshow 1 ) the absence 
of strong Lya emitters (high £W obs ) at high SFR(UV), 2) the 
existence of objects with EW(Lya) b s ^ 100 A for low SFR, 
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3) the possible existence of a maximum value of EW max (Lya) 
as a function of SFR(UV), and 4) a wide range of EW(Lya) b s 
at low star formation rates (SFR(UV) $ 10 M yr _1 ). Although 
the z ~ 3 samples span the largest range in SFR(UV) there does 
not seem to be a systematic trend with redshift. We therefore 
retain all objects, irrespective of their precise z. We suggest the 
following explanations for the trends shown in Fig. [19] 

1) The absence of high EW b s -high SFR objects (1) is due 
to radiation transfer effects and the presence of dust, which un- 
avoidably lead to a stronger reduction of Lya photons compared 
to the adjacent continuum (cf. Sect. |4.7| i. The probability that 
high UV-SFR objects are powered by one or few instantaneous 
bursts is very small; hence their intrinsic Lya- equivalent widths 
must be close to that expected for constant SF, i.e. EW(Lya) D b s > 
70-90 A and up to ~ 250-360 A depending on metallicity (cf. 
Fig. [T31 ). Radiation transfer and the presence of dust reduce 
EW(Lya) b s from the intrinsic to the observed values. 

2) Strong Lya emitters with EW(Lya) b s ^ 100 A corre- 
spond to young 50 Myr) objects, with the exact age limit 
depending on the detailed SF history. Equivalent widths up to 
~ 360 A can be explained by stellar populations with "nor- 
mal" Salpeter-like IMFs and metallicities Z <: 1/50 Z . The 
strong emitters probably show little or no reddening, which 
would otherwise rapidly reduce the observed equivalent width. 
Observationally most of these objects are found as LAE s; few 
LBGs show such strong emission fcf. lShaplev etaT]l2003h . 

3) The trend of EW max (Lya) ("upper envelope") as a func- 
tion of SFR(UV) (i.e. UV magnitude) is most likely due to an 
increase of E(B - V) with SFR(UV). An increasing dust opti- 
cal depth will progressively reduce the intrinsic Lya equivalent 
widths providing a natural continuity between the objects dis- 
cussed in 2) and 1). Presumably the H i column density also 
increases with magnitude (i.e. UV SFR) maintaining a "reason- 
able" spread in the gas-to-dust ratio. Ultimately, the main under- 
lying parameter governing the trends with UV magnitude may 
be the galaxy mass. We will discuss this further below. 

4) The large spread in EW(Lya) Q b s at faint magnitudes 
(SFR(UV) $ 5-10M G yr _1 ) results most likely from two effects: 
first the relatively small amount of dust, which does not elim- 
inate the high EW b s objects, and second the larger variety of 
SF histories/timescales - i.e. an enhanced role of "stochastic SF 
events" - made more plausible for objects of smaller absolute 
scale (mass or total SFR). 

5. Other derived properties 

In other contexts it may be of interest to derive properties such 
as the total neutral hydrogen mass in the outflow and the mass 
outflow rate. The following formulae can be used to estimate 
these quantities. 



5.1. Neutral gas mass 

The H i mass in the shell is related to its column density Nm and 
radius r by 




Assuming r — 1 kpc, we find neutral gas masses of the order of 
~ 2 x 10 6 to 10 s M , with a median of 7 X 10 6 M Q . If the H i were 
found in shells with radii similar to the optical sizes measured 
from D25 for local starbursts, ~ 5-50 kpc, these estimates 



have to increased by a factor 25-2500 . A plau s ible l ower limit 
of r > 1.6/T 1 kpc is derived by Shap lev et al.l (2003) from the 
typical half-light radius of LBGs; similarly the examination of 
close LBG pairs puts an upper limit of r i5 25/z 1 proper kpc on 
the physical dim ension of the gas giving rise to strong interstellar 
absorption lines ( Adelb erger et al.|[2003l) . 



5.2. Mass outflow rates 

Following Petti ni et al. I (120001) the mass loss rate involved in the 
outflow may be estimated as: 




Assuming again r — 1 kpc, we derive mass outflow rates be- 
tween ~ 0.2 and 100 M Q yr . For larger radii Mhi has to be 
increased accordingly. The estimated outflow rates are thus com- 
parable to the star formation rate, as a lready found earlier (e.g. 
Pettin i et al. 2000; [Grimes et al.l f2007). In cases with the lowest 
velocities (quasi-static shells, i.e. in 2-3 of our objects), there is 
most likely no true outflow out of the galaxy. 



6. A unifying scenario for LBGs and LAE, and 
implications 

In paper II, based on our modeling results for cB58 and on em- 
pirical data, we have proposed a unifying scenario to explain 
Lya emission and absorption and the observed trends in LBGs. 
We will now discuss this scenario in light of the results from the 
FDF objects and of the discussion in the present paper, and we 
will show how this scenario should also apply to LAEs. 

In paper II we have suggested that the bulk of the LBGs have 
intrinsically EW(Lya)i nt ~ 60-80 A or larger, and that the main 
physical parameter responsible for the observed variety of Lya 
strengths and profiles in LBGs are A^hi and the accompanying 
variation of the dust content. Here, we propose that the same 
also holds for most LAEs, and that the larger Lya equivalent 
widths found in (some of) these objects are due to younger ages 
of the SF population, quite independently of their SF histories. 
For example, EW(Lya); nt up to ^ 300^-00 A can be expected 
at ages 5 50 Myr and for metallicities ;$ 1/5 Z (cf. Fig. [TBI . 
Any differences between intrinsic and observed Lya equivalent 
widths are then due to radiation transfer and dust effects. 

In paper II, we have shown that these effects transform, e.g. 
for an extinction of E(B-V) ~ 0.3 as found for cB58, a spec- 
trum with an intrinsic Lya emission with EW(Lya)i nt ^ 60 A 
into an absorption-dominated spectrum of LBGs. In the objects 
modeled here, Lya is always found in emission, with observed 
EW(Lya) b s ~ 6-150 A. Again, the intrinsic Lya emission deter- 
mined from our line fits is found to be stronger with EW(Lya); nt 
~ 50-280 A. For all these LBG exa mples, covering the fu ll range 
of Lya line strengths observed bv lShaplev et al.l (120031) . the in- 
trinsic Lya properties are thus compatible with expectations of 
starbursts with constant star formation. For LAEs the same sce- 
nario can also hold, provided the extinction and/or H 1 column 
density are sufficiently low, and provided the age of the highest 
EW(Lya) b s objects is relatively young. 

We now discuss the arguments supporting a continuity and 
even a strong overlap between LBGs and LAE; subsequently we 
will summarise the main observational evidence supporting our 
scenario. 
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Fig. 20. Cumulative distribution of the observed Lya equivalent 
widths of z ~ 3 samples of LAEs and LBGs: Shapley's LBGs 
with R < 25.5 and E W(Lyq)g > 20 A in blue, LAEs from 
iGronwall et ail d2007l) fulfilling the same condition in red, and 
the full LAE sample of Gronwal et al. in black. The probability 
that the blue and red distributions are drawn from the same par- 
ent distribution is 92 %, according to the 2 sided Kolmogorov- 
Smirnov test. 



6.1. LBGs and LAE, overlapping population with a sequence 
driven mostly by mass? 

To be more precise, we suggest the following overlap and dis- 
tinctions between LBGs and LAEs at a given redshift: 

1 . LAEs brighter than a certain limiting magnitude, M^ m , in the 
continuum are the same population as LBGs with strong Lya 
emission (i.e. > EWj'^). 

2. At magnitudes fainter than M\ lm , LAEs represent less mas- 
sive objects than LBGs. 

3. The remaining LBGs, i.e 
Lya in absorption and hence not selected as LAEs, can cover 
a wide range of galaxy masses, SFR, and SF histories 



those with EW(Lya) < EW^ 1 and 



At redshift z ~ 3 the overlap between LAEs and LBGs cor- 
responds to approximately ~ 25% of the LBG population, and 
~ 23 % of the LAE population, as we will discuss below. The 
amount of overlap, hence the value of M]; m varies quite likely 
with redshift. Let us now justify these statements in more detail. 

That there should be a continuity between LBGs and LAEs 
and a strong overlap between the different populations seems 
quite obvious. After all their distinction stems from a "UV- 
continuum" versus "emission line" selection, based, however, 
on the same spectral range (restframe UV). Furthermore the 
Lyman break or drop-out criterion is al so applied as a selec- 
tion criterion for m any LAE samples (e.g. Tanig uchi et alJ2005t 
Ouc hi et al.ll2007l) . and in any case LAEs are expected to show 
a Lyman break, even if too faint to be measured from the cur- 
rent data. A continuity between LBGs and LAEs is also sup- 
ported by the overlap of many of their observed properties , 
such as magnitudes, colors, and others (e.g. lGawiser et al.1l2006t 



ready discussed above (SFR(UV), SFR(Lya), Lya equivalent 
widths). Furthermore, the spatial distribution of LBGs and LAEs 
show also very similar correlation lengths, both at z ~ 3 
and 4.5 where measurements are available dAdelberge r et al.l 
120051: iGawiser et alJl2007t iKovac et al.ll2007h . Finally, at z ~ 3 
the relative number densities of LAE/LBG can also be under- 
stood: in the magnit ude range Ra k < 2 5.5 covered both by 
the L BG sample of IShaplev et alJ d2003l) and the LAE sam- 
ple of IGronwall etal J d2007l) . one has a ratio of LAE/LBG ~ 
1/3, according to the latter authors. This is identical to the frac- 
tion of ~ 25% of LBGs showing strong Lya, EW(Lya) b s ~ 20 
A, since this limit is basically the same as the selection limit 
of the Gronwal samplfl Adopting the same magnitude limit 
(R < 25.5) for both samples, we find also no significant dif- 
ference between the equivalent widths distributions of the z ~ 3 
LBG and LAE sample, as shown in Fig. [20] Indeed, while the 
EW bs distribution of Shapley's total LBG sample differs from 
that of Gronwal, the probability that both distributions are drawn 
from the same parent distribution is 92 % (according to the 
2 sided Kolmogorov-Smirnov test), once this same magnitude 
limit is applied to both samples. This overlap with the LBGs 
concerns 37 out of 160 LAEs (23 %) from Gronwal's statisti- 
cally complete LAE sample at z ~ 3.1. 

We therefore conclude that LAEs above a certain magnitude 
limit represent the same population as LBGs with Lya emission 
(point 1 above); furthermore the LAE selection method allows 
one to find strong emission line objects drawn from a larger 
range of continuum brightness than the LBG selection, as also 
shown in Fig. [19] (point 2). The "remaining" LBGs, correspond- 
ing to three quartiles of Shapley, cannot be found through the 
LAE selection technique, and represent thus a separate class of 
star forming objects. As many LBG studies have shown their 
parameters (mass, SFR etc.) cover a wide range (point 3). 

If a continuity and overlap exists between LBGs and LAE, 
what is then the main parameter(s) "driving" the observed Lya 
strength? From the Lya radiation transfer in expanding shells 
we have seen that the extinction (here described by E(B - V) 
or the dust optical depth respectively) plays a -if not the most 
important- role in determining the strength of Lycfl The next 
important parameter Nui, and the resulting variation in dust/gas 
ratio, as well as possibly deviations from the simple geometry 
considered here, likely introduce the required scatter around an 
oversimplified 1 -dimensional scaling only. Instead of the extinc- 
tion other parameters, such as the stellar mass M* of the galaxy, 
its SFR, or others may be physically more fundamental, as sug- 
gested in paper II. In this case LBGs and LAEs at the same red- 
shift would represent a continuity with decreasing mass, at least 
on average, and correlations between M*, E(B - V), and Nhi 
would naturally explain the observed trends, including those of 
Lya. 

Among LBG samples, iPentericci et alJ d2007l) have e.g. 
found that LBGs at z ~ 4 with Lya emission are less mas- 
sive on average than LBGs with Lya in ab sorption. The same 
trend has also been found by [Erb et al.l d2006l) for z ~ 2 LBGs. 



4 In fact, applying the criteria R AB < 25.5 and EW(Lya)£g > 20 A to 
Shapley's data we find 175 out of 814 LBG, i.e. 21.5 %■ Dropping the 
restric tion on R we obtain 22.5 %, similar to the ~ 25 % of Shaple vet al.1 
J2003h . 

5 Indeed E(B-V) is proportional to the dust absorption optical depth, 
which -associated with the effect of resonant scatterings- determines 
the absorption probability of Lycr photons. 
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Comparing LBGs and LAE, it is well established that LAEs 
show lower UV (i.e. continuum-based) star formation rates, and 
tend to h ave lower stellar masses than LBGs at similar red- 
shift (cf. iReddv et alj|2006t iPirzkal et al|l2007t IGronwall et alj 
l2007tlGawiser et al.ll2007l:IOverzier et al.H2008l:lLai et alj |2008). 
The measured correlation length and galaxy bias of LBGs and 
LAEs are also consistent with a lower mass for the e nsemble of 
LAEs, as discussed recently bv lGawiser et al.l d2007l) . At z ~ 1— 
2 clear correlations between the total SFR and stellar mass M* 
are found for UV star formation rates coveri ng a similar range 
as those observed in z ~ 3 LAEs and LBGs dElbaz et all 2007: 
iDaddi et all 120071: iNoeske et aTll2007l) . Furthermore we know 
that extin ction scales with the total (IR plus UV) SFR (e.g. at 
z ~ 2 see lReddv et alj|2006h . Except for the behaviour of the H i 
column density with galaxy mass, which is so far unknown, all 
quantities show thus the required trend to be consistent with our 
scenario. 

6.2. LBGs and LAEs: age and SF history differences? 

As already pointed out in paper II, stellar populations studies of 
LBGs consistently find ages of several tens to hundreds of Myr 
and favour constant s t ar formation over this timescale (see e.g . 
Ellin gson et all 1 19961: iPettini et al l 120001: Ide Mello et all 120001: 
Papo vich et al.l l200lUShaplev et all 1200 ll) . In this case intrin- 
sic EW(Lya)j nt of ~ 60-100 A are unavoidable (cf. Fig. [15}. 
Although less studies have yet been undertaken on the stel- 
lar populations of LAEs, it appears quite clearly that some of 
them are younger than LBGs. For example, us i ng the sample 
of z = 3.1 LAEs from IGronwall et all d2007l) . iGawiser et ail 
d2007l) : lLai et al.l (120081) find a large spread of ages (from ~ 20 
Myr to ~ 1.6 Gyr) and a decrease of the mean age with decreas- 
ing mass and near-IR luminosity. Similarly, young ages (and low 
mass es) have been found from studies of LAEs a t higher redshift 
(e.g. lPirzkal et al.ll2006tlFinkelstein et al]|2007bl) . 

No age constraint can be derived from Lya alone for objects 
with EW(Lyo') b s 5= 100 A, as such equivalent widths can be 
obtained for constant SF and taking into account various degrees 
of Lya- suppression by dust. On the other hand, EW(Lya) b s > 
100 A requires ages of 55 10-50 Myr at maximum, depending 
on the exact SF history (Fig. [15}; furthermore in such cases the 
amount of dust must be relatively low (e.g. E(B - V) % 0.05- 
0.1 for £W^7EWint = 0.5 and r =0.44-1, according to Eq. 
[4} to maintain a high Lya- equivalent width. In other words, for 
the vast majority of the LBGs Lya provides no age constraint, 
whereas LAEs with EW(Lya) b s ^100 should have ages of $ 
10-50 Myr. These simple predictions appear in good agreement 
with the age estimates of LBGs and LAEs discussed above. 

What can be inferred about the SF timescale/duration in 
LBGs and LAEs and possible differences between the two 
classes? As already mentioned, for the LBGs modeled in de- 
tail here and in paper II, our Lya modeling results require rel- 
atively high intrinsic Lya equivalent widths, corresponding to 
constant star formation over timescales ^ 20-50 Myr or longer; 
Lya alone does not allow us to determine an upper limit for 
the duration of SF. For LAEs with high equivalent widths (^ 
100 A) there is a degeneracy between age and SF duration, 
but long SF timescales are not excluded from their Lya prop- 
erties. Even less strong LAEs could be explained with constant 
SF (and hence unconstrained SF duration). For example objects 
with EW(Lya) bs = 20 A — a value close to the detection cri- 
terion commonly adopted — could be explained with a small 
amount of extinction, E(B - V) gas ~ 0.1, for an intrinsic equiv- 
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Fig. 21. Observed restframe Ly a equivalent w i dth v ersus R AS 
magnitude for the ~ 800 LBGs of lShaplev elaTI (120031) (show as 
crosses). The total EW D b s is plotted here, including both absorp- 
tion and emission components if present. Colors indicate differ- 
ent amounts of extinction E(B - V) estimated from the G - R 
color as described in Shapley et al.: E(B - V) > 0.2 in black, 
red: 0.1-0.2, blue: 0.05-0.1, and yellow: E(B - V) < 0.05. The 
corresponding values of the four LBG groups of Shapley et al. 
are show as filled squares. Note the apparent existence of low 
E(B - V), low EW(Lya) b s objects, which could be dominated 
by stellar absorption. 



alent width three times larger, according to Eq. [4] From Lya 
alone, relatively long SF timescales could thus be possible both 
for the LBGs studied here and for the typical LAEs. However, 
more "bursty", intermittent SF histories may be more plausible 
for lower mass objects, such as LAEs. The suggestion of rela- 
tively long timescales for LAEs is e. g. compatible wit h the ages 
of ~ 160 Myr to 1.6 Gyr derived by lLai et ail (120081) : however 
these age estimates provide probably upper limits, as derived as- 
suming constant SF, and younger ages are also obtained when 
multiple components are allowed for. 

Potentially most interesting to constrain the SF timescale are 
objects with a low Lya equivalent width (faint emission or even 
absorption) and low extinction. If the super-shell model applies 
to such objects and if their velocities are comparable to the typ- 
ical LBGs, observations of Lya in absorption could indicate a 
post-starburst phase dominated by stellar absorption, as clear 
from Fig. 3 of paper II. I n fact, such objects s eem to exist in 
the z ~ 3 LBG sample of Shapl evet al.l 42003), as shown Fig. 
|2T1 Finding objects with shorter SF timescales among the low 
SFR objects, which may also be seen in this Fig., appears more 
plausible given the stronger dependence on "stochasticity" al- 
ready discussed earlier. However, whether these objects are truly 
short duty-cycle, post-starburst objects with little dust remains to 
be verified. Otherwise different geometries and kinematics may 
need to be invoked and the radiation transfer models be reex- 
amined for these objects. Observations of stellar UV features, 
such as Civ and Si iv and other H recombination lines (e.g. Ha) 
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should in principle allow to test this hypothesis and hence pro- 
vide interesting constraints on the duty-cycle of these objects. 

6.3. Open questions and further tests of the scenario 

Although many arguments support our scenario described above, 
several important tests should be carried out. For example, a cor- 
relation between the galaxy mass and the H i column density 
remains to be established. Furthermore, does the dust mass in- 
crease with stellar mass? In other words, how does the dus-to-gas 
ratio evolve with galaxy mass? Also, how well do the proposed 
correlations hold for individual objects? And what governs the 
dispersion around such correlations? 

Concerning LAEs, direct measurements of their ISM prop- 
erties (velocity shifts, ISM covering factors etc.) would be very 
useful. Although it can be expected that these sources show out- 
flows, as well established for the majority of LBGs, such obser- 
vations remain to be done. This could also be important to ex- 
amine possible differences in the homogeneity and the covering 
factor of the ISM in these objects. 

An unclear issue raised by our study concerns the fainter 
LBGs with low Lye equivalent widths or Lya in absorption. Are 
these, presumably less massive objects, sufficiently dust-rich to 
suppress Lya emission, are they post starburst objects, or what 
other process(es) explain their absence of Lya emission? 

A larger number of deep high quality spectroscopic ob- 
servations, multi-wavelength data, and accompanying modeling 
should help to answer these questions in the near future. 

6.4. Implications 

Our scenario proposed for LBGs and LAEs has several notable 
implications, most of them already mentioned in paper II. 

First, our analysis has clearly shown that LBGs exhibit quite 
significant variations in their Lya escape fraction, as apparent 
also from earlier comparisons of Lya and other SFR indicators. 
Such variations must be included in models aimed at understand- 
ing at Lya emission from LAEs and LBGs and alike objects. 

Second, observed distribution functions of EW(Lya) b s and 
the Lya luminosity function must be strongly modified by radia- 
tion transfer and dust effects, and the fraction of objects with low 
EW(Lya) bs (L(Lya)) must be "artificially" overestimated (cf. 
paper II). For a proper interpretation this must be taken into ac- 
count; simple corrections may be applied using e.g. the proposed 
extinction estimate (Eq.[4ji or more sophisticated methods. 

Third, the long/constant star formation rates and "old" ages 
(~ 20-100 Myr) of LBGs derived from stellar population stud- 
ies can be reconciled with the observed Lya properties, when 
radiation transfer and dust are taken into account. In these cir- 
cumstances Lya should basically be independent of the stellar 
population age (measured since the onset of SF); its strength 
and line profile are mostly determined by the extinction and by 
Nui (at least in the case of a homogeneous expanding shell). 
Possibl e trends of Lya emission w i th age, such as indicate d 
e.g. bv IShaplev et alj d200U 120031) ; iFerrara & Ricottil (120061) : 
Pent ericci et alj ( 20071) but apparently found to be conflicting 
among each other, should then be fortuitous. 

Finally, our scenario should naturally explain the in- 
crease of the ratio of LAE/LBG and a higher percentage of 
LBGs with strong Lya emissio n with redshift. These trends 
were suggested earlier (see e.g. iHu et al.1 11998|) and are now 
quite clearly observ ed ( Noll et al]|2004t IShimasaku et al. 2006; 
Nagao et all 120071; lOuchi et all 120071; iDow-Hygelund etall 



120071; iReddv et al.l 120071) . Indeed if extinction decreases with 
increasing z as expected e.g. from the observed metallicity de- 
crease with redshift (cf. [Tremonti et alj 120041 : lErb et al] [2006; 
lAndo et all 120071: iLiu etal.1 12008). our radiation transfer mod- 
els presented here predict stronger Lya emission (see Eq. 0J, 
even if the intrinsic Lya emission properties remain the same. 
Stronger intrinsic emission is additionally expected with de- 
creasing metallicity, due to changes of the ionising/stellar prop- 
erties. Both effects will thus lead to stronger Lya emission and 
hence to an increase of the LAE/LBG ratio with redshift. 



7. Summary and conclusion 

Our 3D Lya and U V continuum radiation transfer code 
dVerhamme et a l. 2006) has been applied to fit the Lya line pro- 
files of 11 LBGs with Lya emission in the FORS Deep Field 
with redshifts between 2. 8 and 5 observed w ith medium spectral 
resolution (R » 2000) by Tap ken et alj d2007l) . The observed line 
profiles show a variety of morphologies, including redshifted 
asymmetric profiles commonly observed, two double-peak pro- 
files, and other more complicated line shapes. 

We have adopted a simple model with a spherically expand- 
ing shell, which consists of neutral hydrogen and mixed in dust, 
surrounding a central starburst emitting a UV continuum plus 
Lya recombination line radiation from its associated H n region. 
With such a model, described by 4 "shell parameters" (the shell 
velocity V exp , its column density Nui, the Doppler parameter b, 
the dust optical depth r a ) and 2 "spectral parameters" (EWj nt and 
FWHMi nt of the input Lya line), we have been able to success- 
fully reproduce all the observed line profiles (see Sect. 01. For 
the majority of the 1 1 LBGs analysed here the medium resolu- 
tion spectra allowed us to constrain these parameters with rel- 
atively few degeneracies. In particular we found that Lya line 
profile fitting alone allows us to determine the amount of dust 
extinction. 

The parameters we have derived for our sample behave as 
follows (see Sect. [4]): 

- The expansion velocity is typically between V exp ~ 150-200 
km s . Higher velocities (300^-00 km s~') may be found 
in two objects. In two out of 1 1 LBGs the neutral medium is 
most likely quasi-static (V e xp ~ 10-25 km s~'); their double- 
peaked line profiles are clearly distinct from the majority of 
the spectra. Our velocity determinations are also supported 
by measurements of interstellar lines in three cases. 

- Except for one object, FDF4691 with a double-peaked Lya 
profile with FWHM ~ 840 km s _1 , all Lya line profiles can 
be reproduced with a relatively low value of the intrinsic Lya 
line width (FWHAf nt ~ 100 km s _1 ); radiation transfer ef- 
fects broaden the line to the observed widths of typically 
FWHM ohs ~ 200-600 km s -1 . We have proposed that the 
observed FWHM(Lya) is related to the H i column density 
to first order. 

- The radial H i column density obtained from our Lya line fits 
ranges from Nm ~ 2 x 10 19 to 7 x 10 20 cirT 2 . Lower values 
may also apply in few objects. We found possible indications 
for an anti-correlation between the observed Lya equivalent 
width and Nm, and for a correlation between FWHM(Lya) 
and Mm- 

- Our Lya profile fits yield values of E(B — V) ~ 0.05-0.2 for 
the gas extinction. These values are also broadly in agree- 
ment with the extinction derived from broad-band SED fits 
for our objects. We found indications for a dust-to-gas ratio 
higher than the Galactic value, and for a substantial scatter. 
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- The escape fraction of Lya photons f e , which follows from 
our Lya profile fits, was found to be determined primarily 
by the extinction. We proposes a simple fit between f e and 
E(B - V) (Eq. O, which for E(B - V) $ 0.3 is also in rea- 
sonable agreement with the Lya escape fraction measured 
in nearby starbursts. Since Lya photons are more strongly 
suppressed by dust than UV continuum photons, a measure- 
ment of EW(Lya) b s can in principle yield E(B - V) (see Eq. 
0), provided the intrinsic Lya equivalent width is known (or 
assumed). 

- Although the LBGs modeled here include objects with ob- 
served Lya equivalent widths down to ~ 6 A in the rest- 
frame (and EW(Lya) obs up to 130-150 A), 10/11 objects 
can clearly be fit with higher intrinsic Lya equivalent widths 
(EW(Lya)i n t ~ 50-100 A or larger), which are expected 
for stellar populations forming constantly over long peri- 
ods (^ 10-100 Myr). In three cases we found indications for 
younger populations. This result agrees with the one found 
in paper II from the modeling of cB58 — an LBG with Lya 
in absorption — and shows that radiation transfer effects and 
dust can lead to a great diversity of observed Lya strengths 
(and profile shapes), even if the intrinsic stellar populations 
and emission line properties of LBGs may be fairly similar. 

- Finally, after correction for radiation transfer and dust effects 
we obtain a better consistency between the use of the Lya 
and UV continuum luminosity as SFR indicators. 

Putting together the model results from this paper (LBG with 
Lya in emission) and those from paper II (modeling of an LBG 
with Lya in absorption) we have been able to reproduce the main 
trends/correlations found in large samples of LBGs and LAEs 
(see Sect. I4.91 . In particular we found that the absence of high 
£Wobs-high SFR objects is due to radiation transfer effects and 
the presence of dust. We also show that the observed correlation 
of velocity shifts betwe en IS lines and Lya with Lya strength 
(cf. IShaplev et al.ll2003l) does not reflect systematic differences 
of the outflow velocities. 

We have proposed that most LBGs and LAEs have intrinsi- 
cally EW(Lya); nt ~ 60-80 A or larger, and that the main physical 
parameter responsible for the observed variety of Lya strengths 
and profiles in LBGs are A/hi and the accompanying variation 
of the dust content. These quantities, in turn, are thought to 
scale mostly with galaxy mass. Such relatively high intrinsic 
equivalent widths are compatible with expectations for rela- 
tively long star formation episodes (k, 10-100 Myr), compatible 
with typical ages and star formation histories derived for LBGs. 
Interestingly we found also indications for relatively low extinc- 
tion (E(B - V) Z> 0. 1 ) LBGs with Lya in absorption in the sample 
of Shapley and collaborators. If true, these could be short duty- 
cycle, post-starburst objects, dominated by stellar Lya absorp- 
tion. 

Analysing the properties of LBGs and LAEs at z ~ 3 we 
have shown that there is a clear overlap between t hese two pop- 
ulation s (Sect. ED: ~ 20-25 % of the LBGs of IShaplev etall 
d2003l) - those with EW^ya)^ > 20 A - correspond to the 
LAEs brighter than Rab = 25.5 mag. The remai ning ~ 77 % of 
the sta tistically complete LAE population from iGronwall et al.l 
(2007) are fainter in the continuum than the LBGs. Radiation 
transfer, dust effects, and changes in the stellar/ionising proper- 
ties should also naturally explain the increase of the LAE/LBG 
ratio, and a higher percentage of LBGs with strong Lya emission 
with increasing redshift. 

Our Lya line profile fitting of LBGs with a detailed radia- 
tion transfer code has shown that we are able to understand a 



diversity of Lya line shapes and strengths and thereby to pro- 
vide interesting quantitative constraints on their gas and stellar 
properties. Most of the observed differences and correlations be- 
tween Lya and other properties can be understood with a simple 
outflow model, where the H i column density and dust content 
(E{B - V)) play the main roles. The same model allow success- 
fully applies to LAE, who appear to be closely related to LBGs. 

Several additional tests of our scenario have been proposed 
and new questions have been raised, e.g. about the star forma- 
tion duration of some LBGs with apparently low extinction (see 
Sect. 16.3b . Numerous interesting implications remain also to be 
worked, e.g. on the interpretation of Lya luminosity functions. 
It is the hope that radiation transfer models, such as those suc- 
cessfully applied here, will help to shed further light on star for- 
mation and galaxy evolution at high redshift. 
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